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Abstract 



In the last few decades the idea came up that by making use of the superposition 
principle from Quantum Mechanics, one can process information in a new and much 
faster way. Hence a new field of information technology, QIT (Quantum Information 
Technology), has emerged. From a physics point of view it is important to find ways of 
implementing these new methods in real systems. One of the most basic tasks required 
for QIT is the ability to connect different components of a Quantum Computer by 
quantum wires that obey the superposition principle. Since superpositions can be 
very sensitive to noise this turns out to be already quite difficult. Recently, it was 
suggested to use chains of permanently coupled spin-1/2 particles (quantum chains) 
for this purpose. They have the advantage that no external control along the wire is 
required during the transport of information, which makes it possible to isolate the 
wire from sources of noise. The purpose of this thesis is to develop and investigate 
advanced schemes for using quantum chains as wires. We first give an introduction to 
basic quantum state transfer and review existing advanced schemes by other authors. 
We then introduce two new methods which were created as a part of this thesis. First, 
we show how the fidelity of transfer can be made perfect by performing measurements 
at the receiving end of the chain. Then we introduce a scheme which is based on 
performing unitary operations at the end of the chain. We generalise both methods 
and discuss them from the more fundamental point of view of mixing properties of a 
quantum channel. Finally, we study the effects of a non-Markovian environment on 
quantum state transfer. 
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Notation 



X,Y,Z Pauli matrices 

X n ,Y n , Z n Pauli matrices acting on the Hilbert-space of qubit n 

|0), |1) Single qubit state in the canonical basis 

|0) Quantum chain in the product state |0) <8> ■ • • <8> |0) 

|n) "Single excitation" state X n \0) 

Trx Partial trace over subsystem X 

1 1 ... 1 1 Euclidean vector norm 

1 1 ... 1 1 1 Trace norm 

1 1 ... 1 1 2 Euclidean matrix norm 



We also use the following graphical representation: 

|o> = 
|i> = ® 

= a\0) + P\l)= O 



o)-0000000 
n)-OOCXs)000 

nth qubit 




sender controlled region: 
( Alice ) quantum gates 

and measurements coupling 



uncontrolled region receiver 
("Bob") 
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1 Introduction 



The Hilbert space that contains the states of quantum mechanical objects is huge, 
scaling exponentially with the number of particles described. In 1982, Richard Feyn- 
man suggested to make use of this as a resource for simulating quantum mechanics in a 
quantum computer, i.e. a device where the physical interaction could be "programmed" 
to yield a specific Hamiltonian. This has led to the new fields of Quantum Computa- 
tion and Quantum Information. A quantum computer can solve questions one could 
never imagine to solve using an ordinary computer. For example, it can factorise 
large numbers into primes efficiently, a task of greatest importance for cryptography. 
It may thus be a surprise that more than twenty years after the initial ideas, these 
devices still haven't been built or only in ridiculously small size. The largest quantum 
computer so far can only solve problems that any child could solve within seconds. 
A closer look reveals that the main problem in the realisation of quantum comput- 
ers is the "programming", i.e. the design of a specific (time-dependent) Hamiltonian, 
usually described as a set of discrete unitary gates. This turns out to be extremely 
difficult because we need to connect microscopic objects (those behaving quantum 
mechanically) with macroscopic devices that control the microscopic behaviour. Even 
if one manages to find a link between the micro- and the macroscopic world, such as 
laser pulses and electric or magnetic fields, then the connection introduces not only 
control but also noise (dissipation and decoherence) to the microscopic system, and 
its quantum behaviour is diminished. 

The vision of this thesis is to develop theoretical methods narrowing the gap between 
what is imagined theoretically and what can be done experimentally. As a method 
we consider chains (or more general graphs) of permanently coupled quantum systems. 
This idea has been originally put forward by S. Bose for the specific task of quantum 
communication [1]. Due to the permanent coupling, these devices can in principle be 
built in such a way that they don't require external control to perform their tasks, 
just like a mechanical clockwork. This also overcomes the problem of decoherence as 
they can be separated from any source of noise. Unfortunately, most schemes that 
have been developed so far still require external control, though much less than an 
"ordinary" quantum computer. Furthermore, internal dispersion in these devices is 
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leading to a decrease of their fidelity. A third problem is, that for building these 
devices the permanent couplings still need to be realised, although only once, and 
experimental constraints such as resolution and errors need to be considered. We are 
thus left with the following questions: which is the best way to perform quantum 
state transfer using a permanently coupled graph? How much control do we need, 
and how difficult will it be to implement the couplings? How do errors and noise 
affect the scheme? All these points are highly related and it cannot be expected to 
find an absolute, i.e. system independent answer. The purpose of this research is to 
develop advanced schemes for the transfer of quantum information, to improve and 
generalise existing ideas, to relate them to each other and to investigate their stability 
and efficiency. 



1.1 Quantum Computation and Quantum Information 

In this Section we review some of the basic concepts of Quantum Computation. We 
will be very brief and only focus on those aspects that we require later on in the thesis. 
A more detailed introduction can be found in [2]. 

In information science, an algorithm is a list of instructions that a computer performs 
on a given input to achieve a specific task. For instance, a factoring algorithm has an 
arbitrary integer as its input, and gives its prime factors as an output. A quantum 
factoring algorithm can be thought of in a similar way, i.e. it has an integer as input, 
and its prime factors as an output. In-between however it encodes information in a 
quantum mechanical system. Due to the superposition principle, the information of a 
quantum system cannot be represented as bits. The valid generalisation of the bit to 
the quantum case is called qubit. The possible states of a qubit are written as 

a\0)+P\l), (1.1) 

where a, (3 are normalised complex coefficients, and |0) and |1) are vectors of a two- 
dimensional complex vector space. Peter W. Shor has shown in a famous paper [3] 
that the detour of representing the intermediate part of a factoring algorithm in a 
quantum system (as well as using quantum gates, see below) can be very beneficial: 
it runs much faster. This is important, because many cryptographic methods rely 
on factoring algorithms being slow. Shor's algorithm is definitely not the only reason 
why it would be very nice to have a quantum computer, i.e. a machine that represents 
information in a quantum way and can perform instructions on it, and many more 
details can be found in the textbook mentioned above. 
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Algorithms on a computer can be represented as list of logical operations on bits. 
Likewise, a (standard) quantum algorithm can be represented as a list of quantum 
logical operations, or quantum gates, acting on qubits. The most general quantum 
algorithm is given by an arbitrary unitary operator. A universal set of gates is a 
set such that any quantum algorithm (i.e. unitary operator) can be decomposed 
into a sequence of gates belonging to this set. In the standard model of quantum 
computation, one assumes that such a set is available on the machine [4]. Also the 
ability to perform measurements is assumed. We refer to this as the full control case. 

From a information theoretic point of view, qubits are not only useful objects to 
perform algorithms with, but also very interesting from a fundamental point of view. 
To give a (too simple) analogy consider the following. If you read the word "chocolate", 
you can associate a positive/negative or neutral feeling of whether you would like to eat 
some chocolate now. However, what was the state of your mind concerning chocolate 
before you read the word? Unless you were already craving for chocolate beforehand, 
or you have just eaten a lot, your mind was probably undecided. Moreover, it would 
have been very difficult - if not impossible - to describe to someone in plain language 
which opinion you had about the chocolate before you read the word. 

In a similar manner, the quantum information contained in a single arbitrary and 
unknown qubit cannot be described by classical information. When it is measured, it 
behaves like a normal bit in the sense that the outcome is only or 1, but when it 
is not measured, it behaves in some way as if it was undecided between and 1. Of 
course one has to be very careful with these analogies. But for the purpose of this 
thesis it is important to stress that quantum information cannot be transported by 
any classical methods [5]. This is why it is so important and also so difficult to develop 
new wires, dubbed quantum wires, that are capable of doing this. 

1.2 Quantum state transfer along short distances 

In theory, additional devices for the transfer of unknown quantum states are not 
required for building a quantum computer, unless it is being used for typical quantum 
communication purposes, such as secret key distribution [4]. This is because the 
universal set of gates on the quantum computer can be used to transfer quantum 
states by applying sequences of two-qubit swap gates (Fig. 11.1( 1. 

However in practice it is crucial to minimise the required number of quantum gates, 
as each gate typically introduces errors. In this light it appears costly to perform A — 1 
swap gates between nearest neighbours to just move a qubit state over a distance of 
A sites. For example, Shor's algorithm on A qubits can be implemented by only 
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Figure 1.1: In areas of universal control, quantum states can easily be transferred by 
sequences of unitary swap gates Sj t k between nearest neighbours. 

log N quantum gating operations [6] if long distant qubit gates are available. These 
long distant gates could consist of local gates followed by a quantum state transfer. 
If however the quantum state transfer is implemented as a sequence of local gates, 
then the number of operations blows up to the order of N gates. The quantum state 
transfer can even be thought of as the source of the power of quantum computation, as 
any quantum circuit with log N gates and local gates only can be efficiently simulated 
on a classical computer [7,8]. 

A second reason to consider devices for quantum state transfer is related to scala- 
bility. While small quantum computers have already been built [9], it is very difficult 
to build large arrays of fully controllable qubits. A black box that transports unknown 
quantum states could be used to build larger quantum computers out of small compo- 
nents by connecting them. Likewise, quantum state transfer can be used to connect 
different components of a quantum computer, such as the processor and the memory 
(see also Fig. 11.21) . On larger distances, flying qubits such as photons, ballistic electrons 
and guided atoms/ions are considered for this purpose [10,11]. However, converting 
back and forth between stationary qubits and mobile carriers of quantum information 
and interfacing between different physical implementations of qubits is very difficult 
and worthwhile only for short communication distances. This is the typical situation 
one has to face in solid state systems, where quantum information is usually contained 
in the states of fixed objects such as quantum dots or Josephson junctions. In this case 
permanently coupled quantum chains have recently been proposed as prototypes of re- 
liable quantum communication lines [1,12]. A quantum chain (also referred to as spin 
chain) is a one-dimensional array of qubits which are coupled by some Hamiltonian 
(cf. Fig. II. 3p . These couplings can transfer states without external classical control. 
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In many cases, such permanent couplings are easy to build in solid state devices (in 
fact a lot of effort usually goes into suppressing them). The qubits can be of the same 
type as the other qubits in the device, so no interfacing is required. 



Input 



Controller 



Quantum 
processor 



Quantum 
memory 



Output 



Figure 1.2: Schematic layout of a quantum computer. The solid arrows represent the 
flow of quantum information, and the dashed arrows the flow of classical information. 




Figure 1.3: Permanently coupled quantum chains can transfer quantum states without 
control along the line. Note that the ends still need to be controllable to initialise and 
read out quantum states. 

Another related motivation to consider quantum chains is that they can simplify the 
layout of quantum devices on wafers. A typical chip can contain millions of qubits, and 
the fabrication of many qubits is in principle no more difficult than the fabrication of 
a single one. In the last couple of years, remarkable progress was made in experiments 
with quantum dots [13,14] and super-conducting qubits [15,16]. It should however 
be emphasised that for initialisation, control and readout, those qubits have to be 
connected to the macroscopic world (see Fig- 1 1 .2j) . For example, in a typical flux qubit 
gate, microwave pulses are applied onto specific qubits of the sample. This requires 
many (classical) wires on the chip, which is thus a compound of quantum and classical 
components. The macroscopic size of the classical control is likely to be the bottleneck 
of the scalability as a whole. In this situation, quantum chains are useful in order to 
keep some distance between the controlled quantum parts. A possible layout for such 
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a quantum computer is shown in Fig. 11.41 It is built out of blocks of qubits, some 
of which are dedicated to communication and therefore connected to another block 
through a quantum chain. Within each block, arbitrary unitary operations can be 
performed in a fast and reliable way (they may be decomposed into single and two- 
qubit operations). Such blocks do not currently exist, but they are the focus of much 
work in solid state quantum computer architecture. The distance between the blocks 
is determined by the length of the quantum chains between them. It should be large 
enough to allow for classical control wiring of each block, but short enough so that 
the time-scale of the quantum chain communication is well below the time-scale of 
decoherence in the system. 



\ooo/ \oocv* 

'OOOOOOOOOOOOOOOOOOOO' 

/ooo\ /ooo\ 



Figure 1.4: Small blocks (grey) of qubits (white circles) connected by quantum chains. 
Each block consists of (say) 13 qubits, 4 of which are connected to outgoing quantum 
chains (the thick black lines denote their nearest-neighbour couplings). The blocks are 
connected to the macroscopic world through classical wires (thin black lines with black 
circles at their ends) through which arbitrary unitary operations can be triggered on 
the block qubits. The quantum chains require no external control. 

Finally, an important reason to study quantum state transfer in quantum chains 
stems from a more fundamental point of view. Such systems in principle allow tests of 
Bell-inequalities and non-locality in solid-state experiments well before the realisation 
of a quantum computer. Although quantum transport is quite an established field, 
the quantum information point of view offers many new perspectives. Here, one looks 
at the transport of information rather than excitations, and at entanglement [17-20] 
rather than correlation functions. It has recently been shown that this sheds new 
light on well-known physical phenomena such as quantum phase transitions [21-24], 
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quantum chaos [25-28] and localisation [29,30]. Furthermore, quantum information 
takes on a more active attitude. The correlations of the system are not just calculated, 
but one also looks at how they may be changed. 

1.3 Implementations and experiments 

As we have seen above, the main advantage of state transfer with quantum chains is 
that the qubits can be of the same type as those used for the quantum computation. 
Therefore, most systems that are thought of as possible realisations of a quantum 
computer can also be used to build quantum chains. Of course there has to be some 
coupling between the qubits. This is typically easy to achieve in solid state sys- 
tems, such as Josephson junctions with charge qubits [31,32], flux qubits [33,34] (see 
also Fig. I1.5p or quantum dots dots using the electrons [35,36] or excitons [37,38]. 
Other systems where quantum chain Hamiltonians can at least be simulated are NMR 
qubits [39-41] and optical lattices [42]. Such a simulation is particularly useful in the 
latter case, where local control is extremely difficult. Finally, qubits in cavities [43,44] 
and coupled arrays of cavities were considered [45,46]. 




Figure 1.5: A quantum chain consisting of N = 20 flux qubits [34] (picture and exper- 
iment by Floor Paauw, TU Delft). The chain is connected to four larger SQUIDS for 
readout and gating. 

For the more fundamental questions, such as studies of entanglement transfer, non- 
locality and coherent transport, the quantum chains could also be realised by systems 
which are not typically thought of as qubits, but which are natural spin chains. These 
can be molecular systems [47] or quasi-lD solid state materials [48,49]. 

1.4 Basic communication protocol 

We now review the most basic transport protocol for quantum state transfer, initially 
suggested in [1]. For the sake of simplicity, we concentrate on the linear chain setting, 
though more general graphs of qubits can be considered in the same way. The protocol 
consists of the following steps: 
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1. Initialise the quantum chain in the ground state 

\G). (1.2) 

2. Put an arbitrary and unknown qubit with (possibly mixed) state p at the sending 
end of the chain 

P ^Tn{\G)(G\}. (1.3) 

3. Let the system evolve under its Hamiltonian H for a time t 

exp {-iHt} p ® Tri { | G) (G\ } exp {iHt} . (1.4) 

4. Pick up the quantum state at the end of the chain 

c = Tri i ... iA r_i [exp{-iHt}p®Tri{\G){G\}exp{iHt}] . (1.5) 

Some practical aspects how to realise these steps are discussed in the next section. 
For the moment, we will concentrate on the quality of quantum state transfer given 
that the above steps can be performed. From a quantum information perspective, the 
above equations describe a quantum channel [5] r that maps input states p at one end 
of the chain to output states r(p) = a on the other end. A very simple measure of the 
quality of such a quantum channel is the fidelity [2,50,51] 

F(p,a) = (rtsjp^ap 1 / 2 ) . (1.6) 

More advanced measures of the quality of transfer will be discussed in Chapter [3l Note 
also that some authors define the fidelity without taking the square of the trace. It is 
a real- valued, symmetric function with range between and 1, assuming unity if and 
only if p = a. Since the transported state that is an unknown result of some quantum 
computation, we are interested in the minimal fidelity 

F = mm p F(p,r(p)). (1.7) 

We remark that some authors also assume an equal distribution of input states and 
compute the average fidelity [1]. Using the strong concavity of the fidelity [2] and the 
linearity of r we find that the minimum must be assumed on pure input states, 

F = min^|T(^M. (1-8) 
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In the present context, Fq = Fo(H,t) is a function of of the Hamiltonian H of the 
quantum chain (through the specific role of the ground state in the protocol and 
through the time evolution), and of the time interval t that the system is evolving in 
the third step of the protocol. 

1.4.1 Initialisation and end-gates 

There are two strong assumptions in the protocol from the last section. The first one 
is that the chain can be initialised in the ground state \G). How can that be achieved 
if there is no local control along the chain? The answer appears to be quite easy: one 
just applies a strong global magnetic field and strong cooling (such as laser cooling or 
dilution refrigeration) and lets the system reach its ground state by relaxation. The 
cooling needs to be done for the remaining parts of the quantum computer anyway, 
so no extra devices are required. However there is a problem with the time-scale of 
the relaxation. If the system is brought to the ground state by cooling, it must be 
coupled to some environment. But during the quantum computation, one clearly does 
not want such an environment. This is usually solved by having the time-scale of the 
computation much smaller (say microseconds) than the time-scale of the cooling (say 
seconds or minutes). But if the quantum chain should be used multiple times during 
one computation, then how is it reset between each usage? This is important to avoid 
memory effects [52], and there are two solutions to this problem. Either the protocol 
is such that at the end the chain is automatically in the ground state. Such a protocol 
usually corresponds to perfect state transfer. The other way is to use the control at 
the ends of the chain to bring it back to the ground state. A simple cooling protocol 
is given by the following: one measures the state of the last qubit of the chain. If it 
is in |0), then one just lets the chain evolve again and repeats. If however it is found 
to be in |1), one applies the Pauli operator X to flip it before evolving and repeating. 
It will become clear later on in the thesis that such a protocol typically converges 
exponentially fast to the ground state of the chain. 

The second assumption in the last section is that the sender and receiver are capable 
of swapping in and out the state much quicker than the time-scale of the interaction of 
the chain. Alternatively, it is assumed that they can switch on and off the interaction 
between the chain and their memory in such time-scale. It has recently been shown [33] 
that this is not a fundamental problem, and that finite switching times can even slightly 
improve the fidelity if they are carefully included in the protocol. But this requires to 
solve the full time-dependent Schrodinger equation, and introduces further parameters 
to the model (i.e. the raise and fall time of the couplings). For the sake of simplicity, 
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we will therefore assume that the end gates are much faster then the time evolution 
of the chain (see also Section 

1.4.2 Symmetries 

The dimensionality of the Hilbert space TL of a quantum chain of N qubits is 2 N . 
This makes it quite hopeless in general to determine the minimal fidelity Eq. (|1.8p for 
long quantum chains. Most investigations on quantum state transfer with quantum 
chains up to date are therefore concentrating on Hamiltonians with additional sym- 
metries. With few exceptions [21,22,34,53] Hamiltonians that conserve the number of 
excitations are considered. In this case the Hilbert space is a direct sum of subspaces 
invariant under the time evolution, 

N 

n = @H t , (1.9) 

with dimWf = (A, and where t is the number of excitations. These Hamiltonians are 
much easier to handle both analytically and numerically, and it is also easier to get 
an intuition of the dynamics. Furthermore, they occur quite naturally as a coupling 
between qubits in the relevant systems. We stress though that there is no fundamental 
reason to restrict quantum chain communication to this case. 

1.4.3 Transfer functions 

The space TCo only contains the state |0) which is thus always an eigenstate of H. We 
will assume here that it is also the ground state, 

\G) = |0>. (1.10) 

This can be achieved by applying a strong global magnetic field (or equivalent) to the 
system. The space Hi is spanned by the vectors {\k), k = 1, . . . ,N} having exactly 
one excitation. The above protocol becomes: 

1. Initialise the quantum chain in the ground state 

|0> (1.11) 

2. Put an arbitrary and unknown qubit in the pure state \ip) = a\0) + at the 
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sending end of the chain 

a\0)+P\l) (1.12) 



3. Let the system evolve for a time t 

a\0) + (3exp{-iHt} |1) (1.13) 



4. Pick up the quantum state at the end of the chain (see [1]) 

rW = (1 -p(t))|0)<0| +p(tM<# (1.14) 

with the minimal fidelity given by 

F = min^(V|r(V)|V) (1-15) 
= p(t) + (1 - p(t))mm f |(0|V)| 2 = p(t). (1.16) 



The function p(t) is the transition probability from the state |1) to \N) given by 

p(t) = \(N\exp{-iHt}\l)\ 2 . (1.17) 

We see that in the context of quantum state transfer, a single parameter suffices to 
characterise the properties of an excitation conserving chain. The averaged fidelity [1] 
is also easily computed as 

Jpjtj p(t) 1 
F= + 6 + 2' (L18) 

Even more complex measures of transfer such as the quantum capacity only depend on 
pit) [54]. It is also a physically intuitive quantity, namely a particular matrix element 
of the time evolution operator, 

/n,m(t) = H exp {-iHt} \m) (1.19) 

= Y,e- iE k t (n\E k )(E k \m), (1.20) 
k 

where \Ek) and E k are the eigenstates and energy levels of the Hamiltonian in 7i\. 
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1.4.4 Heisenberg Hamiltonian 



The Hamiltonian chosen in [1] is a Heisenberg Hamiltonian 

N-l N 

'J 



H = ~!zJ2 ( X n x n+i + Y n Y n+l + Z n Z n+1 ) -BY,Z n + c, (1.21) 



n=l 



n=l 



with a constant term 



J(N - 1) 



+ NB 



(1.22) 



added to set the ground state energy to 0. For J > it fulfils all the assumptions 
discussed above, namely its ground state is given by |0) and it conserves the number 
of excitations in the chain. The Heisenberg interaction is very common and serves 
here as a typical and analytically solvable model for quantum state transfer. 



In the first excitation subspace H\, the Heisenberg Hamiltonian Eq. (|1.21j) is ex- 
pressed in the basis {\n}} as 



/ 1 



-1 



\ 



V 



1 1 



(1.23) 



/ 



A more general study of such tridiagonal matrices can be found in a series of articles 
on coherent dynamics [55-58]. Some interesting analytically solvable models have also 
been identified [56,57,59] (we shall come back to that point later). 



For the present case, the eigenstates of Eq (I1.23P are [1] 



1 + ^fcO 
N 



N 

£ 

n=l 



COS 



irk 
2N' 



(2n-l) 



\n) (k = 0, 



N-l) 



(1.24) 



with the corresponding energies given by 

E k = 2B + 2 J 



1-cos- 



(1.25) 



The parameter B has no relevance for the fidelity but determines the stability of the 
ground state (the energy of the first excited state is given by 2B). The minimal fidelity 
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for a Heisenberg chain is given by 



p(t) = N~ 



1 + E exp {-2, Jt(l - cos ^) J (-l) fc (l + cos ^ 



As an example, Fig 11.61 shows p(t) for N = 50. 

0.14 




(1.26) 



_i i i i i y i v 

" 10 20 30 40 50 60 70 80 90 100 

Time [1/J] 

Figure 1.6: Minimal fidelity p(t) for a Heisenberg chain of length N = 50. 



1.4.5 Dynamic and Dispersion 

Already in [1] has been realised that the fidelity for quantum state transfer along spin 
chains will in general not be perfect. The reason for the imperfect transfer is the 
dispersion [60] of the information along the chain. Initially the quantum information 
is localised at the sender, but as it travels through the chain it also spreads (see Fig. II . 71 
and Fig. II. 8p . This is not limited to the Heisenberg coupling considered here, but a 
very common quantum effect. Due to the dispersion, the probability amplitude peak 
that reaches Bob is typically small, and becomes even smaller as the chains get longer. 



The fidelity given Eq. (I1.26P is shown in Fig. 11.61 We can see that a wave of 
quantum information is travelling across the chain. It reaches the other end at a time 
of approximately 

*peak ~ ( L2? ) 
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Figure 1.7: Snapshots of the time evolution of a Heisenberg chain with N = 50. Shown 
is the distribution |/ n ,i(t)| 2 of the wave-function in space at different times if initially 
localised at the first qubit. 
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Figure 1.8: Mean and variance of the state |1) as a function of time. Shown is the case 
N = 50 with the y-axis giving the value relative to the mean N/2 + 1 and variance 
(iV 2 - 1)/12 of an equal distribution -4= ^ \n). 
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As a rough estimate of the scaling of the fidelity with respect to the chain length 
around this peak we can use [1,61] (see also Fig. 1 1 . 9[) 

| /jVil (i)| 2 « |2J^(|)| 2 « | (^) V3 ai[(!) (N - §)]| 2 , (1.28) 

where Jjv(x) is a Bessel function of first kind and ai{x) is the Airy function. The airy 
function ai{x) has a maximum of 0.54 at x = —1.02. Hence we have 

K*peak) = \fN,i(^j)\ 2 « L82iV- 2 / 3 . (1.29) 

It is however possible to find times where the fidelity of the chain is much higher. 
The reason for this is that the wave-packet is reflected at the ends of the chain and 
starts interfering with itself (Fig ll.6jl . As the time goes on, the probability distribution 
becomes more and more random. Sometimes high peaks at the receiving end occur. 
From a theoretical point of view, it is interesting to determine the maximal peak 
occurring, i.e. 

Pm(T) = max pit). (1.30) 

0<t<T 

As we can see in Fig. 11.101 there is quite a potential to improve from the estimate 
Eq. (fl~29l) . 




1 - 1 — 1 1 1 

20 22 24 26 28 30 

Time [1/J] 

Figure 1.9: Approximation of the transfer amplitude for N = 50 around the first 
maximum by Bessel and Airy functions [1,61]. 
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Figure 1.10: pm(T) as a function of T for different chain lengths. The solid curve is 
given by 1.82(2T)~ 2//3 and corresponds to the first peak of the probability amplitude 
(Eq. ESP 

We will now show a perhaps surprising connection of the function pm(T) to number 
theory. Some speculations on the dependence of the fidelity on the chain length being 
divisible by 3 were already made in [1], but not rigorously studied. As it turns out, for 
chains with prime number length the maximum of the fidelity is actually converging 
to unity (see Fig. II. 10) . To show this, we first prove the following 

Lemma 1.1 Let N be an odd prime. Then the set 

jcos— (fc = 0,l,...„— ^— )| (1.31) 

is linear independent over the rationals Q. 
Proof Assume that 

N-l 

^A fc cos^ = (1.32) 

k=0 

with Afc G Q. It follows that 

N-l 

Ef ikir —ikn 1 

k 1 exp ~/v~ + exp ~ /v - I = ^' 33 ^ 

fc=o ^ ' 
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and hence 



JV-l 
2 



JV-l 
2 



i(N - k)ir 



A fc exp — - 2^ A fe exp - = 0. 

fc=0 fc=o 



N 



Changing indexes on the second sum we get 



jV— 1 
2 



iV 



Afc exp — - 2^ *N-k exp = 0. 



fc=0 



. JV+1 



N 



(1.34) 



(1.35) 



and finally 



where 



Afc exp — = 0, 

fc=0 



(1.36) 



Ao — 2Ao 

Afc = Afc (k = 1, 



N- 1, 



' 2 ' 
iV + 1 

Afc = -A^_fc (fc = _2^,...,jv-l). 



(1.37) 
(1.38) 

(1.39) 



Since N is prime, the roots of unity in Eq. (I1.36P are all primitive and therefore linearly 
independent over Q [62, Theorem 3.1, p. 313]. Hence Afc = for all k. ■ 



Theorem 1.1 (Half 


recurrence) Let N be an odd prime. 


For a Heisenberg 


chain of length N we 


have 










lim Pm(T) = lim 


max p(t) 
0<t<T 


= 1. 


(1.40) 



Proof The eigenfrequencies of the Hamiltonian in the first excitation sector Hi are 
given by 



E k = 2B + 2 J 



Trk 

1-cos- 



0,1, 



N -11 



(1.411 



Using Kronecker's theorem [63] and Lemma [TTT1 the equalities 



exp{itE k } = (-1) e 



k 2(B+J)t 



(A: = 0,1,...,, 



iV- 1, 



(1.42) 
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can be fulfilled arbitrarily well by choosing an appropriate t. Since 

cos— = -cos V - ' , (1.43) 

the equalities (|1.42|) are then also fulfilled arbitrarily well for k = 0, . . . , N — 1. This 
is known as as sufficient condition for perfect state transfer in mirror symmetric 
chains [64], where the eigenstates can be chosen such that they are alternately sym- 
metric and antisymmetric. Roughly speaking, Eq. (|1.42|) introduces the correct phases 
(a sign change for the antisymmetric eigenstates) to move the state |1) to \N) and 
hence the theorem. ■ 

Remark 1.1 The time-scale for finding high valued peaks is however exponential in 
the chain length [63]. Therefore the above theorem has little practical use. For non- 
prime chain lengths, the eigenfrequencies are not sufficiently independent to guarantee 
a perfect state transfer, with the algebraic dimensionality of the roots of unity for non- 
prime N given by the Euler totient function <f){N) [62, Theorem 3.1, p. 313]. We also 
remark that due to its asymptotic character, the above result is not contradicting [65], 
where it was shown that chains longer than N > 4 never have perfect fidelity. 

Having proved that there are many chains that can in principle perform arbitrarily well, 
it is important to find a cut-off time for the optimisation Eq. (|1.30|) . Faster transfer 
than linear in N using local Hamiltonians is impossible due to the Lieb-Robinson 
bound [66,67], which is a "speed limit" in non-relativistic quantum mechanics giving 
rise to a well defined group velocity. Transport faster than this group velocity is 
exponentially suppressed. Going back to the motivation of quantum state transfer, a 
natural comparison [37] for the time-scale of quantum state transfer is given by the 
time it would take to perform a sequence of swap gates (cf . Fig 11.1)1 that are realised 
by a pairwise switchable coupling Hamiltonian 

^(X n X n+l +Y n Y n+l ). (1.44) 

This time is linear in the chain length: 

Wp = (1.45) 

Ideally one could say that the time for quantum state transfer should not take much 
longer than this. However one may argue that there is a trade-off between quick 
transfer on one hand, and minimising control on the other hand. A second cut-off 
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time may be given by the decoherence time of the specific implementation. But short 
decoherence times could always be counteracted by increasing the chain coupling J. 
A more general and implementation independent limit is given by the requirement 
that the peak width At^ ea ^ should not be too small with respect to the total time. 
Otherwise it is difficult to pick up the state at the correct time. For the first peak, 
we can estimate the width by using the full width at half height of the airy function. 
From Eq. (|1.28p we get an absolute peak width of Atp ea ^ 0.72A rl / 3 /J and a relative 
width of 

« 1.44AT 2 /3. (1 46) 

^peak 

This is already quite demanding from an experimental perspective and we conclude 
that the transfer time should not be chosen much longer than those of the first peak. 



1.4.6 How high should p(t) be? 

We have not discussed yet what the actual value of p(t) should be to make such a 
spin chain useful as a device for quantum state transfer. p(t) = corresponds to no 
state transfer, p(t) = 1 to a perfect state transfer. But what are the relevant scales 
for intermediate p(t)? In practice, the quantum transfer will suffer from additional 
external noise (Chapter [7J and also the quantum computer itself is likely to be very 
noisy. From this point of view, requiring p(t) = 1 seems a bit too demanding. 

From a theoretical perspective, it is interesting that for any p(t) > 0, one can al- 
ready do things which are impossible using classical channels, namely entanglement 
transfer and distillation [2]. The entanglement of formation between the sender (^4/- 
ice) and the receiver (Bob) is simply given by \/p(t) [1]. This entanglement can be 
partially distilled [68] into singlets, which could be used for state transfer using tele- 
portation [2]. It is however not known how much, i.e. at which rate, entanglement 
can be distilled (we will develop lower bounds for the entanglement of distillation in 
Section [231 and Section f3~4^ . Also, entanglement distillation is a quite complex pro- 
cedure that requires local unitary operations and measurements, additional classical 
communication, and multiple chain usages; and few explicit protocols are known. This 
is likely to preponderate the benefits of using a quantum chain. 

When the chain is used without encoding and further operation, the averaged fidelity 
Eq. (|1.18p becomes better than the classical averaged fidelity [1] when p(t) > 3 — 

x By "classical fidelity", we mean the fidelity that can be achieved by optimising the following pro- 
tocol: Alice performs measurements on her state and sends Bob the outcome through a classical 
communication line. Bob then tries to rebuild the state that Alice had before the measurement 
based on the information she sent. For qubits, the classical fidelity is given by 2/3 [69]. 
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2y2. Following the conclusion from the last subsection that the first peak is the most 
relevant one, this would mean that only chains with length until N = 33 perform 
better than the classical fidelity. 

Finally, the quantum capacity [54, 70] of the channel becomes non-zero only when 
p(t) > 1/2, corresponding to chain lengths up to N = 6. Roughly speaking, it is a 
measure of the number of perfectly transmitted qubits per chain usage that can be 
achieved asymptotically using encoding and decoding operations on multiple channel 
usages. The quantum capacity considered here is not assumed to be assisted by a 
classical communication, and the threshold of p(t) > 0.5 to have a non-zero quantum 
capacity is a result of the non-cloning theorem [2]. This is not contradicting the 
fact that entanglement distillation is possible for any p(t) > 0, as the entanglement 
distillation protocols require additional classical communication. 

All the above points are summarised in Fig. 11.111 We can see that only very short 
chains reach reasonable values (say > 0.6) for the minimal fidelity. 
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Figure 1.11: Quantum capacity, entanglement of formation (EOF), a lower bound for 
the entanglement of distillation (EOD) and the averaged fidelity as a function of p(t). 
We also show the corresponding chain length which reaches this value as a first peak 
and the classical threshold 3 — 2\/2- The explicit expression for the quantum capacity 
plotted here is given in [54], and the lower bound of the entanglement of distillation 
will be derived in Section 13.41 
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1.5 Advanced communication protocols 

We have seen in the last section that without much further effort, i.e. entanglement 
distillation, unmodulated Heisenberg chains are useful only when they are very short. 
Shortly after the initial proposal [1] it has been shown that there are ways to achieve 
even perfect state transfer with arbitrarily long chains. These advanced proposals can 
roughly be grouped into four categories, which we will now briefly describe. 

1.5.1 Engineered Hamiltonians 

The Heisenberg model chosen by Bose features many typical aspects of coherent trans- 
port, i.e. the wave-like behaviour, the dispersion, and the almost-periodicity of the 
fidelity. These features do not depend so much on the specific choices of the parame- 
ters of the chain, such as the couplings strengths. There are however specific couplings 
for quantum chains that show a quite different time evolution, and it was suggested 
in [71] and independently in [72] to use these to achieve a perfect state transfer: 



These values for engineered couplings also appear in a different context in [57,73]. The 
time evolution under the Hamiltonian (|1.47p features an additional mirror symmetry: 
the wave-packet disperses initially, but the dispersion is reversed after its centre has 
passed the middle of the chain (Fig. II. 121) , This approach has been extended by various 
authors [19,53,59,64,65,74-84], and many other choices of parameters for perfect or 
near perfect state transfer in various settings were found [59,81,83]. 

1.5.2 Weakly coupled sender and receiver 

A different approach of tuning the Hamiltonian was suggested in [85]. There, only 
the first and the last couplings j of the chain are engineered to be much weaker than 
the remaining couplings J of the chain, which can be quite arbitrary. The fidelity 
can be made arbitrarily high by making the edge coupling strengths smaller. It was 
shown [86, 87] that to achieve a fidelity of 1 — 5 in a chain of odd length, it takes 
approximately a time of 



JV-l 




(1.47) 



n=l 



2Ntt/V6 



(1.48) 
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Figure 1.12: Snapshots of the time evolution of a quantum chain with engineered 
couplings (I1.47P for N = 50. Shown is the distribution of the wave-function in space 
at different times if initially localised at the first qubit (compare Fig. II. 7|) . 
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and the coupling ratio has to be approximately j/J « yS/N . Some specific types of 
quantum chains which show high fidelity for similar reasons were also investigated [88- 
91], 

1.5.3 Encoding 

We have seen in Subsec. 11.4.61 that if p(t) < 1/2, the fidelity cannot be improved 
by using any encoding/decoding strategy (because the quantum capacity is zero). 
However it is possible to change the protocol described in Sec. 11.41 slightly such that 
the fidelity is much higher. This can be thought of as a "hardware encoding", and 
was suggested first in [60]. There, it was assumed that the chain consists of three 
sections: one part of length 2./V 1 / 3 controlled by the sending party, one "free" part 
of length N and one part of length « 2.8./V 1 / 3 controlled by the receiving party. The 
sender encodes the qubit not only in a single qubit of the chain, but in a Gaussian- 
modulated superposition of his qubits. These Gaussian packets are known to have 
minimal dispersion. Likewise, the receiver performs a decoding operation on all qubits 
he controls. Near-perfect fidelity can be reached. 

1.5.4 Time-dependent control 

Finally, a number of authors found ways of improving the fidelity by time-dependent 
control of some parameters of the Hamiltonian. In [92] it is shown that if the end 
couplings can be controlled as arbitrary (in general complex valued) smooth functions 
of time the encoding scheme [60] could be simulated without the requirement of ad- 
ditional operations and qubits. Another possibility to achieve perfect state transfer 
is to have an Ising interaction with additionally pulsed global rotations [40,93,94]. 
Further related methods of manipulating the transfer by global fields were reported 
in [25,28,95-98]. 

1.6 Motivation and outline of this work 

While the advanced transfer protocols have shown that in principle high fidelity can 
be achieved with arbitrarily long chains, they have all come at a cost. Engineering 
each coupling of the Hamiltonian puts extra demands on the experimental realisation, 
which is often already at its very limits just to ensure the coherence of the system. 
Furthermore, the more a scheme relies on particular properties of the Hamiltonian, the 
more it will be affected by imperfections in its implementation [84,99]. For example, 
simulating an engineered chain of length N = 50 with a (relative) disorder of 5%, we 
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get a fidelity peak of 0.95 ± 0.02. For a disorder of 10% we get 0.85 ± 0.05. The weakly 
coupled system is very stable for off-site disorder [85], but suffers strongly from on-site 
disorder (i.e. magnetic fields in z— direction) at the ends of the chain. For example, 
for a chain of N = 50 with edge couplings j = 0.01 and the remaining couplings 
being J = 1, we find that already a magnetic field of the order of 0.00001 lowers the 
fidelity to 0.87 ± 0.12. For fields of the order of 0.00005 we find 0.45 ± 0.32. This is 
because these fluctuations must be small with respect to the small coupling, so there 
is a double scaling. Also, the time-scale of the transfer is longer than in other schemes 
(note though that this may sometimes even be useful for having enough time to pick up 
the received state). On the other hand, encoding and time-dependent control require 
additional resources and gating operations. It is not possible to judge independently of 
the realisation which of the above schemes is the "most practical" one. We summarise 
the different aspects that are important in the following five criteria for quantum state 
transfer: 

1. High efficiency: How does the fidelity depend on the length of the chain? Which 
rate [74,81,100] can be achieved? 

2. Minimal control: How many operations are required to achieve a certain fidelity? 
WhereQ is control required? 

3. Minimal resources: What additional resources are required? 

4. Minimal design: How general is the coupling typJH? What values of the coupling 
strengths are allowed? 

5. Robustness: How is the fidelity affected by static disorder, by time-dependent 
disorder, by gate and timing errors, and by external noise such as decoherence 
and dissipation? 

At the start of this research, only the engineering and encoding schemes were available. 
The engineering schemes are strong in the points 2 and 3, but quite weak in the points 
4 and 5. The encoding scheme on the other hand has its weakness in points 2 and 3. 
It was hence desirable to develop more balanced schemes. Since most experiments 
in Quantum Information are extremely sensitive and at the cutting edge of their 
parameters (i.e. requiring extremely low temperatures, well tuned lasers, and so forth, 

2 For example, gates at the ends of the chain are always needed for write-in and read-out, and thus 
"cheaper" than gates along the chain. Global control along the whole chain is often easier than 
local control. 

3 Often the coupling type is already fixed by the experiment 
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to maintain their quantum behaviour), we particularly wanted to find schemes which 
are strong in the points 4 and 5. Also, from a more fundamental point of view, we 
were interested in seeing how much information on the state of a quantum chain could 
be obtained by the receiver in principle, and how the receiver might even be able to 
prepare states on the whole chain. 

The main achievements of this thesis are two schemes for the transfer of quantum 
information using measurements (Chapter [2] and [3]) or unitary operations (Chapter [5] 
andE]) at the receiving end of the chain. Since both schemes use convergence properties 
of quantum operations, it seemed natural to investigate these properties in a more 
abstract way (Chapter [4]). There, we found a new way of characterising mixing maps, 
which has applications beyond quantum state transfer, and may well be relevant for 
other fields such as chaos theory or statistical physics. Finally, in Chapter [7] we 
discuss problems quantum state transfer in the presence of external noise. The results 
in Chapters [3]|6] were developed in collaboration with Vittorio Giovannetti from Scuola 
Normale, Pisa. Much of the material discussed in this thesis has been published or 
submitted for publication [101-110]. 
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2.1 Introduction 

The role of measurement in quantum information theory has become more active 
recently. Measurements are not only useful to obtain information about some state 
or for preparation, but also, instead of gates, for quantum computation [111]. In the 
context of quantum state transfer, it seems first that measurements would spoil the 
coherence and destroy the state. The first indication that measurements can actually 
be used to transfer quantum information along anti-ferromagnetic chains was given 
in [24]. However there the measurements had to be performed along the whole chain. 
This may in some cases be easier than to perform swap gates, but still requires high 
local accessibility. We take a "hybrid" approach here: along the chain, we let the 
system evolve coherently, but at the receiving end, we try to help the transfer by 
measuring. The main disadvantage of the encoding used in the protocols above is 
that once the information dispersed, there is no way of finding out where it is without 
destroying it. A dual rail encoding [112] as used in quantum optics on the other 
hand allows us to perform parity type measurements that do not spoil the coherence 
of the state that is sent. The outcome of the measurement tells us if the state has 
arrived at the end (corresponding to a perfect state transfer) or not. We call this 
conclusively perfect state transfer. Moreover, by performing repetitive measurements, 
the probability of success can be made arbitrarily close to unity. As an example of 
such an amplitude delaying channel, we show how two parallel Heisenberg spin chains 
can be used as quantum wires. Perfect state transfer with a probability of failure lower 
than P in a Heisenberg chain of N qubits can be achieved in a time-scale of the order 
of 0.33J _1 A^ L7 | lnP|. We demonstrate that our scheme is more robust to decoherence 
and non-optimal timing than any scheme using single spin chains. 

We then generalise the dual rail encoding to disordered quantum chains. The scheme 
performs well for both spatially correlated and uncorrelated fluctuations if they are 
relatively weak (say 5%). Furthermore, we show that given a quite arbitrary pair of 
quantum chains, one can check whether it is capable of perfect transfer by only local 
operations at the ends of the chains, and the system in the middle being a black box. 
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We argue that unless some specific symmetries are present in the system, it will be 
capable of perfect transfer when used with dual rail encoding. Therefore our scheme 
puts minimal demand not only on the control of the chains when using them, but also 
on the design when building them. 

This Chapter is organised as follows. In Section f2T2l we suggest a scheme for quan- 
tum communication using two parallel spin chains of the most natural type (namely 
those with constant couplings). We require modest encodings (or gates) and measure- 
ments only at the ends of the chains. The state transfer is conclusive, which means 
that it is possible to tell by the outcome of a quantum measurement, without destroy- 
ing the state, if the transfer took place or not. If it did, then the transfer was perfect. 
The transmission time for conclusive transfer is not longer than for single spin chains. 
In Section 12.31 we demonstrate that our scheme offers even more: if the transfer was 
not successful, then we can wait for some time and just repeat the measurement, 
without having to resend the state. By performing sufficiently many measurements, 
the probability for perfect transfer approaches unity. Hence the transfer is arbitrarily 
perfect. We will show in Section 12.41 that the time needed to transfer a state with a 
given probability scales in a reasonable way with the length of the chain. In Section 
12.51 we show that encoding to parallel chains and the conclusiveness also makes our 
protocol more robust to decoherence (a hitherto unaddressed issue in the field of quan- 
tum communication through spin chains). In the last part of this chapter, we show 
how this scheme can be generalised to disordered chains (Sections l2.6ti2.lQ]) and even 
coupled chains (Section l2.1ip . 

2.2 Scheme for conclusive transfer 

We intend to propose our scheme in a system-independent way with occasional refer- 
ences to systems where conditions required by our scheme are achieved. We assume 
that our system consists of two identical uncoupled spin-l/2-chains (1) and (2) of 
length N, described by the Hamiltonian 

H = ® 1^ + /(D ® - Egl® I {2 \ (2.1) 

The term identical states that and are the same apart from the label of the 
Hilbert space they act on. The requirement of parallel chains instead of just one is 
not a real problem, since in many experimental realisations of spin chains, it is much 
easier to produce a whole bunch of parallel uncoupled [48,49] chains than just a single 
one. 
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A B 

Figure 2.1: Two quantum chains interconnecting A and B. Control of the systems is 
only possible at the two qubits of either end. 

We assume that the ground state of each chain is \0) i} i.e. a ferromagnetic ground 
state, with H^' \0) i = E g \0) i , and that the subspace consisting of the single spin 
excitations \n) i is invariant under Let us assume that the state that Alice wants 
to send is at the first qubit of the first chain, i.e. 

|^> 1 = a|0) 1 + /3|l) 1) (2.2) 

and that the second chain is in the ground state 1 0)2 • The aim of our protocol is to 
transfer quantum information from the 1st ("Alice") to the Ath ("Bob") qubit of the 
first chain: 

\iI> a ) x ^\M>b)i = <*\Q) 1 + P\N) 1 . (2.3) 

The first step (see also Fig. I2.2p is to encode the input qubit in a dual rail [112] by 
applying a NOT gate on the first qubit of system (2) controlled by the first qubit of 
system (1) being zero, resulting in a superposition of excitations in both systems, 

|a(0))=a|0,l) + /9|l,0), (2.4) 

where we have introduced the short notation |n, m) = \n)\ (g> \rn)2- This is assumed 
to take place in a much shorter time-scale than the system dynamics. Even though a 
2-qubit gate in solid state systems is difficult, such a gate for charge qubits has been 
reported [15]. For the same qubits, Josephson arrays have been proposed as single spin 
chains for quantum communication [31]. For this system, both requisites of our scheme 
are thus available. In fact, the demand that Alice and Bob can do measurements and 
apply gates to their local qubits (i.e. the ends of the chains) will be naturally fulfilled 
in practice since we are suggesting a scheme to transfer information between quantum 
computers (as described in Section [L2l) . 
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|o) (2) — ef 



spin chain (1) 
spin chain (2) 




wait again if 



Figure 2.2: Quantum circuit representation of conclusive and arbitrarily perfect state 
transfer. The first gate at Alice's qubits represents a NOT gate applied to the second 
qubit controlled by the first qubit b eing zero. The qubit 1^^)^ on the left hand side 
represents an arbitrary input state at Alice's site, and the qubit \iPb)i represents the 
same state, successfully transferred to Bob's site. The i^-gate represents the unitary 
evolution of the spin chains for a time interval of t£. 



Under the system Hamiltonian, the excitation in Eq. (|2.4|) will travel along the two 
systems. The state after the time t\ can be written as 



N 



\<j>(h)) = Y J fnAti)\s(n)}, (2.5) 



n=l 



where \s(n)) = a |0, n) + (3\n, 0) and the complex amplitudes f n ,i(ti) are given by 
Eq. (|1.19|) . We can decode the qubit by applying a CNOT gate at Bob's site. Assuming 
that this happens on a time-scale much shorter than the evolution of the chain, the 
resulting state is given by 



N-l 



^2 fn,i(h) \s(n)} + f^ih) W B ) 1 |7V> 2 ■ (2-6) 



n=l 



Bob can now perform a measurement on his qubit of system (2). If the outcome of 
this measurement is 1, he can conclude that the state has been successfully 

transferred to him. This happens with the probability \fN,i(ti)\ 2 . If the outcome is 
0, the system is in the state 

TV— 1 

= £ /n,i(ti) \s(n)) , (2.7) 



n=l 



where -P(l) = 1 — |/Ar,i(ti)| 2 is the probability of failure for the first measurement. If 
the protocol stopped here, and Bob would just assume his state as the transferred one, 
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the channel could be described as an amplitude damping channel [54], with exactly the 
same fidelity as the single chain scheme discussed in [1]. Note that here the encoding is 
symmetric with respect to a and (3, so the minimal fidelity is the same as the averaged 
one. 

But success probability is more valuable than fidelity: Bob has gained knowledge 
about his state, and may reject it and ask Alice to retransmit (this is known as a 
quantum erasure channel [113]). Of course in general the state that Alice sends is the 
unknown result of some quantum computation and cannot be sent again easily. This 
can be overcome in the following way: Alice sends one e-bit on the dual rail first. If 
Bob measures a success, he tells Alice, and they both start to teleport the unknown 
state. If he measures a failure, they reset the chains and start again. Since the joint 
probability of failure converges exponentially fast to zero this is quite efficient. In fact 
the conclusive transfer of entanglement is possible even on a single chain by using the 
same chain again instead of a second one [114]. This can be seen as a very simple 
entanglement distillation procedure, achieving a rate of \fN,i(t)\ 2 /2. However the chain 
needs to be reset between each transmission (see Section 11.4.11 for problems related 
to this), and Alice and Bob require classical communication. We will show in the 
next section, that the reuse of the chain(s) is not necessary, as arbitrarily perfect state 
transfer can already achieved in the first transmission. 



2.3 Arbitrarily perfect state transfer 

Because Bob's measurement has not revealed anything about the input state (the 
success probability is independent of the input state), the information is still residing 
in the chain. By letting the state (|2.7I) evolve for another time t2 and applying the 
CNOT gate again, Bob has another chance of receiving the input state. The state 
before performing the second measurement is easily seen to be 

1 - 

-j== {/*,l(*2 + *l) " fn,N(t2)fN,l(ti)} \s(n)} . (2.8) 
V ^K L ) n=l 

Hence the probability to receive the qubit at Bobs site at the second measurement is 

|/jv,i(^ + ti)-/jv,7v(t 2 )/7V,i(ti)| 2 . (2.9) 



P(l) 

If the transfer was still unsuccessful, this strategy can be repeated over and over. 
Each time Bob has a probability of failed state transfer that can be obtained from the 
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generalisation of Eq. (|2.8|l to an arbitrary number of iterations. The joint probability 
that Bob fails to receive the state all the time is just the product of these probabilities. 
We denote the joint probability of failure for having done I unsuccessful measurements 
as P(£). This probability depends on the time intervals tg between the [l — l)th and 
£th measurement, and we are interested in the case where the tg are chosen such that 
the transfer is fast. It is possible to write a simple algorithm that computes P(£) 
for any transition amplitude f r ,s(t)- Figure E31 shows some results for the Heisenberg 
Hamiltonian given by Eq. (|1.2ip . 
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Figure 2.3: Semilogarithmic plot of the joint probability of failure P{t) as a function of 
the number of measurements £. Shown are Heisenberg spin-l/2-chains with different 
lengths N. The times between measurements tg have been optimised numerically. 



An interesting question is whether the joint probability of failure can be made 
arbitrarily small with a large number of measurements. In fact, the times tg can be 
chosen such that the transfer becomes arbitrarily perfect. We will prove this in the 
next Chapter, where a generalisation of the dual rail scheme and a much wider class 
of Hamiltonians is considered. In the limit of large number of measurements, the spin 
channel will not damp the initial amplitude, but only delay it. 
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2 « 1.82iV~ 2/3 . (2.11) 



2.4 Estimation of the time-scale the transfer 

The achievable fidelity is an important, but not the only criterion of a state transfer 
protocol. In this Section, we give an heuristic approach to estimate the time that it 
needs to achieve a certain fidelity in a Heisenberg spin chain. The comparison with 
numeric examples is confirming this approach. 

Let us first describe the dynamic of the chain in a very qualitative way. Once Alice 
has initialised the system, an excitation wave packet will travel along the chain. As 
shown in Subsection 11.4.51 it will reach Bob at a time of the order of 

N 

*peak ~ 2j> ( 2,1 °) 

with an amplitude of 

/•V,l(*peak) 

It is then reflected and travels back and forth along the chain. Since the wave packet 
is also dispersing, it starts interfering with its tail, and after a couple of reflections 
the dynamic is becoming quite randomly. This effect becomes even stronger due to 
Bobs measurements, which change the dynamics by projecting away parts of the wave 
packet. We now assume that 2£p ea ^ (the time it takes for a wave packet to travel 
twice along the chain) remains a good estimate of the time-scale in which significant 
probability amplitude peaks at Bobs site occur, and that Eq. (|2.1ip remains a good 
estimate of the amplitude of these peak£|. Therefore, the joint probability of failure 
is expected to scale as 

P{1) « (l - 1.82iV" 2/3 ) (2.12) 

in a time of the order of 

t(l) « 2t max £ = J- X NL (2.13) 

If we combine Eq. (|2.12p and (|2.13p and solve for the time t(P) needed to reach a 
certain probability of failure P, we get for N 1 

t(P) « 0.55 J- 1 N 5/3 |ln F| . (2.14) 

We compare this rough estimate with exact numerical results in Fig. 12.41 The best fit 



lr This is not a strong assumption. If the excitation was fully randomly distributed, the probability 
would scale as iV -1 . By searching for good arrival times, this can be slightly increased to N~ 2 ^ 3 . 
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for the range shown in the figure is given by 

~~' " (2-15) 



t(P) = 0.33J _1 iV 5/3 |lnP| 



We can conclude that the transmission time for arbitrarily perfect transfer is scaling 
not much worse with the length N of the chains than the single spin chain schemes. 
Despite of the logarithmic dependence on P, the time it takes to achieve high fidelity 
is still reasonable. For example, a system with N = 100 and J = 20K * ks will take 
approximately 1.3ns to achieve a fidelity of 99%. In many systems, decoherence is 
completely negligible within this time-scale. For example, some Josephson junction 
systems [115] have a decoherence time of 500ns, while trapped ions have even 

larger decoherence times. 




Figure 2.4: Time t needed to transfer a state with a given joint probability of failure 
P across a chain of length N. The points denote exact numerical data, and the fit is 
given by Eq. (|2.15p . 
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2.5 Decoherence and imperfections 

If the coupling between the spins J is very small, or the chains are very long, the 
transmission time may no longer be negligible with respect to the decoherence time. 
It is interesting to note that the dual rail encoding then offers some significant general 
advantages over single chain schemes. Since we are suggesting a system-independent 
scheme, we will not study the effects of specific environments on our protocol, but just 
qualitatively point out its general advantages. 

At least theoretically, it is always possible to cool the system down or to apply a 
strong magnetic field so that the environment is not causing further excitations. For 
example in flux qubit systems, the system is cooled to « 2hmK to ensure that the 
energy splittingA 3> ksT [116]. Then, there are two remaining types of quantum 
noise that will occur: phase noise and amplitude damping. Phase noise is a serious 
problem and arises here only when an environment can distinguish between spin flips 
on the first chain and spin flips on the second chain. It is therefore important that the 
environment cannot resolve their difference. In this case, the environment will only 
couple with the total z-component 



of the spins of both chains at each position n. This has been discussed for spin-boson 
models in [117,118] but also holds for spin environments as long as the chains are close 
enough. The qubit is encoded in a decoherence-free subspace [119] and the scheme is 
fully robust to phase noise. Even though this may not be true for all implementations 
of dual rail encoding, it is worthwhile noticing it because such an opportunity does 
not exist at all for single chain schemes, where the coherence between two states with 
different total z-component of the spin has to be preserved. Having shown one way of 
avoiding phase noise, at least in some systems, we now proceed to amplitude damping. 

The evolution of the system in presence of amplitude damping of a rate V can be 
easily derived using a quantum-jump approach [120]. This is based on a quantum 
master equation approach, which is valid in the Born-Markov approximation [121] 
(i.e. it holds for weakly coupled environments without memory effects). Similarly to 
phase noise, it is necessary that the environment acts symmetrically on the chains. 
The dynamics is then given by an effective non-Hermitian Hamiltonian 



z (i) + z m 



(2.16) 




(2.17) 



n 
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if no jump occurs. If a jump occurs, the system is back in the ground state |0). The 
state of the system before the first measurement conditioned on no jump is given by 

N 

e- rt Y,fn,i(t)\s(n)), (2.18) 

n=l 

and this happens with the probability of e~ 2rt (the norm of the above state). If a 
jump occurs, the system will be in the ground state 

\/l ~ e~ 2n |0, 0} . (2.19) 

The density matrix at the time t is given by a mixture of (|2.18|) and (|2.19|) . In case 
of (|2.19p . the quantum information is completely lost and Bob will always measure 
an unsuccessful state transfer. If Bob however measures a success, it is clear that no 
jump has occurred and he has the perfectly transferred state. Therefore the protocol 
remains conclusive, but the success probability is lowered by e~ 2Tt . This result is 
still valid for multiple measurements, which leave the state (|2.19p unaltered. The 
probability of a successful transfer at each particular measurement £ will decrease 
by e~ 2Tt ^\ where t{£) is the time at which the measurement takes place. After a 
certain number of measurements, the joint probability of failure will no longer decrease. 
Thus the transfer will no longer be arbitrarily perfect, but can still reach a very high 
fidelity. Some numerical examples of the minimal joint probability of failure that can 
be achieved, 

oo 

lim P(£) m TT ( 1 - 1.357V- 2/3 e-^) (2.20) 

1=1 

are given in Fig. 12.51 For J/T = 50.K" ns nearly perfect transfer is still possible for 
chains up to a length of N 40. 

Even if the amplitude damping is not symmetric, its effect is weaker than in single 
spin schemes. This is because it can be split in a symmetric and asymmetric part. 
The symmetric part can be overcome with the above strategies. For example, if the 
amplitude damping on the chains is Ti and T2 with T\ > T2, the state (|2.18p will be 

N 

/n.l(*) { ae ' T2t 1°' n > + P e ~ Tlt l n ' °>} ( 2 - 21 ) 

n=l 

N 

« e- T2t Y,fn,i{t)\s{n)) (2.22) 

n=l 
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provided that t < (Ti - T 2 ) 1 . Using a chain of length N = 20 with J = 20K * k B 
and T^ 1 = 4ns, T^ 1 = 4.2ns we would have to fulfil t <C 164ns. We could perform 
approximately 10 measurements (cf. Eq. (|2.13|) ) without deviating too much from the 
state (|2.22|) . In this time, we can use our protocol in the normal way. The resulting 
success probability given by the finite version of Eq. (|2.20|) would be 75%. A similar 
reasoning is valid for phase noise, where the environment can be split into common 
and separate parts. If the chains are close, the common part will dominate and the 
separate parts can be neglected for short times. 
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Figure 2.5: The minimal joint probability of failure P(i) for chains with length N in 
the presence of amplitude damping. The parameter J/T of the curves is the coupling 
of the chain (in Kelvin) divided by the decay rate (ns _1 ). 



2.6 Disordered chains 

The main requirement for perfect transfer with dual rail encoding in the above is 
that two identical quantum chains have to be designed. While this is not so much 
a theoretical problem, for possible experimental realizations of the scheme [31] the 
question arises naturally how to cope with slight asymmetries of the channels. We 
are now going to demonstrate that in many cases, perfect state transfer with dual rail 
encoding is possible for quantum chains with differing Hamiltonians. 

By doing so, we also offer a solution to another and perhaps more general problem: 
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if one implements any of the schemes for quantum state transfer, the Hamiltonians will 
always be different from the theoretical ones by some random perturbation. This will 
lead to a decrease of fidelity in particular where specific energy levels were assumed 
(see [84, 99] for an analysis of fluctuations affecting the engineered chains described in 
Subsection II. 5. This problem can be avoided using the scheme described below. In 
general, disorder can lead to a Anderson localisation [29,30,122] of the eigenstates (and 
therefore to low fidelity transport of quantum information). In this section however 
this is not relevant, as we consider only short chains (N < 100) and small disorder 
(« 10% of the coupling strength), and the localisation length is much longer then 
the length of the chain. We will show numerically that the dual rail scheme can still 
achieve arbitrarily perfect transfer for a uniformly coupled Heisenberg Hamiltonian 
with disordered coupling strengths (both for the case of spatially correlated and un- 
corrected disorder). Moreover, for any two quantum chains, we show that Bob and 
Alice can check whether their system is capable of dual rail transfer without directly 
measuring their Hamiltonians or local properties of the system along the chains but 
by only measuring their part of the system. 



2.7 Conclusive transfer in the presence of disorder 

We consider two uncoupled quantum chains (1) and (2), as shown in Fig. 12.61 The 

chains are described by the two Hamiltonians and with total Hamiltonian 
given by 

H = H m ®I {2) ® H {2) , (2.23) 
and the time evolution operator factorising as 

U(t) = exp (-itf (1) t) ® exp (-iH^t) . (2.24) 

For the moment, we assume that both chains have equal length N, but it will become 
clear in Section I2U1 that this is not a requirement of our scheme. All other assumptions 
remain as in the first part of the chapter. 

Initially, Alice encodes the state as 

a\0, 1) +(3\1,0). (2.25) 
This is a superposition of an excitation in the first qubit of the first chain and an 
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A 



B 



Figure 2.6: Two disordered quantum chains interconnecting A and B. Control of the 
systems is only possible at the two qubits of either end. 



excitation in the first qubit of the second chain. The state will evolve into 

N 

Y,{a9n,i{t)\0,n) + (3f njl {t)\n,0)} , (2.26) 

71=1 

with 

f n>1 (t) = (n,0\U(t)\l,0) (2.27) 

g nA (t) = (0,n|tf(t)|0,l>. (2.28) 

In Section l2~2l these functions were identical. For differing chains this is no longer the 
case. We may, however, find a time t\ such that the modulus of their amplitudes at 
the last spins are the same (see Fig. 12.71) . 

9N,i(h) = e^f NA (t 1 ). (2.29) 



At this time, the state (I2.26P can be written as 

JV-l 

{«ffn,i(*i) |0, n) + /3/ n> i(ti) \n, 0)} + 

n=l 

fN,i(h) {e^a |0, N) + | AT, 0)} . (2.30) 
Bob decodes the state by applying a CNOT gate on his two qubits, with the first qubit 
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Figure 2.7: The absolute values of the transition amplitudes /jv,i(i) and gN,i(t) for 
two Heisenberg chains of length N = 10. The couplings strengths of both chains were 
chosen randomly from the interval [0.8J, 1.2J] . The circles show times where Bob can 
perform measurements without gaining information on a and (3. 



as the control bit. The state thereafter is 
JV-l 



{«ffn,i(ti) |0, n) + /3/ n ,i(ti) \n, 0)} + 

n=l 

f N ,l(h) {e^a |0) (1) + |AT)«} ® |iV> (2) . 



(2.31) 



Bob then measures his second qubit. Depending on the outcome of this measurement, 
the systems will either be in the state 



N-l 



— V {a<7 n ,i(ii) |0, n) + Pf n ,x{t x ) \n, 0)} 



or m 



{e^a|0} (1) +/3|JV) (1) }®|iV) (2) , 



(2.32) 



(2.33) 



where p\ = 1— |/Ar,i(ti)| 2 = 1 — |5jv,i(*i)| 2 is the probability that Bob has not received 
the state. The state (|2.33l) corresponds to the correctly transferred state with a known 
phase error (which can be corrected by Bob using a simple phase gate). If Bob finds the 
system in the state (I2.32p . the transfer has been unsuccessful, but the information is 
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still in the chain. We thus see that conclusive transfer is still possible with randomly 
coupled chains as long as the requirement (|2.29l) is met. This requirement will be 
further discussed and generalised in the next section. 



2.8 Arbitrarily perfect transfer in the presence of disorder 

If the transfer was unsuccessful, the state (|2.32p will evolve further, offering Bob 
further opportunities to receive Alice's message. For identical quantum chains, leads 
to a success for any reasonable Hamiltonian (Section l3.6p . For differing chains, this is 
not necessarily the case, because measurements are only allowed at times where the 
probability amplitude at the end of the chains is equal, and there may be systems where 
this is never the case. In this section, we will develop a criterion that generalises Eq. 
(I2.29P and allows to check numerically whether a given system is capable of arbitrarily 
perfect state transfer. 

The quantity of interest for conclusive state transfer is the joint probability P{£) 
that after having checked I times, Bob still has not received the proper state at his end 
of the chains. Optimally, this should approach zero if £ tends to infinity. In order to 
derive an expression for P(£), let us assume that the transfer has been unsuccessful for 
£ — 1 times with time intervals te between the the £th and the (£ — l)th measurement, 
and calculate the probability of failure at the £th measurement. In a similar manner, 
we assume that all the £ — 1 measurements have met the requirement of conclusive 
transfer (that is, Bob's measurements are unbiased with respect to a and f3) and derive 
the requirement for the ith measurement. 

To calculate the probability of failure for the £t\i measurement, we need to take 
into account that Bob's measurements disturb the unitary dynamics of the chain. If 
the state before a measurement with the outcome "failure" is , the state after the 
measurement will be 




(2.34) 



where Q is the projector 

Q = I - |0, N) (0, N\ - | AT, 0) (N, 0| , (2.35) 

and pi is the probability of failure at the Ith measurement. The dynamics of the 
chain is alternating between unitary and projective, such that the state before the ith 
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measurement is given by 



l e 

j== II {U(tk)Q} {a\l,Q)+/3 10, 1)} 



where 



P(t-l) = l[pk. 



(2.36) 



(2.37) 



Note that the operators in (I2.36P do not commute and that the time ordering of the 
product (the index k increases from right to left) is important. The probability that 
there is an excitation at the iVth site of either chain is given by 



P(£ 



a\ 2 \F(£)\ 2 + \3\ 2 



\G(£)\ 2 } 



with 



and 



F(£) = (N,0\l[{U(t k )Q}\l,0) 



k=l 

e 



G(£) = (0,N\l[{U(t k )Q}\0, 1) 



k=l 



(2.38) 



(2.39) 



(2.40) 



Bob's measurements are therefore unbiased with respect to a and f3 if and only if 



|F(*)| = |G(*)| VI 



(2.41) 



In this case, the state can still be transferred conclusively (up to a known phase). The 
probability of failure at the £th measurement is given by 



Pt = l 



\F(£)\ 2 
P{£-1) 



(2.42) 



It is easy (but not very enlightening) to show [103] that the condition (I2.4ip is equiv- 
alent to 



Y[{U(t k )Q} \1,0) 



k=l 



k=l 



VI, 



(2.43) 



and that the joint probability of failure - if at each measurement the above condition 
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is fulfilled - is simply given by 



2 



P{£)= \[{U{t k )Q}\l,Q) 



(2.44) 



k=l 



It may look as if Eq. (|2.43p was a complicated multi-time condition for the measuring 
times ti, that becomes increasingly difficult to fulfil with a growing number of mea- 
surements. This is not the case. If proper measuring times have been found for the 
first i — 1 measurements, a trivial time that fulfils Eq. (|2.43l) is = 0. In this 
case, Bob measures immediately after the (£ — l)th measurement and the probability 
amplitudes on his ends of the chains will be equal - and zero (a useless measurement). 
But since the left and right hand side of Eq. (|2.43|) when seen as functions of tg are 
both almost-periodic functions with initial value zero, it is likely that they intersect 
many times, unless the system has some specific symmetry or the systems are com- 
pletely different. Note that we do not claim at this point that any pair of chains will 
be capable of arbitrary perfect transfer. We will discuss in the next system how one 
can check this for a given system by performing some simple experimental tests. 

2.9 Tomography 

Suppose someone gives you two different experimentally designed spin chains. It may 
seem from the above that knowledge of the full Hamiltonian of both chains is necessary 
to check how well the system can be used for state transfer. This would be a very 
difficult task, because we would need access to all the spins along the channel to 
measure all the parameters of the Hamiltonian. In fact by expanding the projectors in 
Eq. (|2.43p one can easily see that the only matrix elements of the evolution operator 
which are relevant for conclusive transfer are 



Physically, this means that the only relevant properties of the system are the transition 
amplitudes to arrive at Bob's ends and to stay there. The modulus of /jv,i(i) and 
fN,N(t) can be measured by initialising the system in the states |1, 0) and \N, 0) and 
then performing a reduced density matrix tomography at Bob's site at different times 



fN,l(t) 
/iV,Ar(i) 
9N,l(t) 

9N,N(t) 



(N,0\U(t)\l,0) 
(N,0\U{t)\N, 0) 
(0,JV|Z7(t)|0,l) 
(0,N\U(t)\0,N) 



(2.45) 
(2.46) 
(2.47) 
(2.48) 
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t, and the complex phase of these functions is obtained by initialising the system 
in (|0, 0) + |1, 0)) /y/2 and (|0, 0) + \N, 0)) /y/2 instead. In the same way, g N ,i(t) 
and gN,N(t) are obtained. All this can be done in the spirit of minimal control at 
the sending and receiving ends of the chain only, and needs to be done only once. 
It is interesting to note that the dynamics in the middle part of the chain is not 
relevant at all. It is a black box (see Fig. 12. 8|) that may involve even completely 
different interactions, number of spins, etc., as long as the total number of excitations 
is conserved. Once the transition amplitudes [Equations (|2.45|) - (|2.48p ] are known, one 





Figure 2.8: The relevant properties for conclusive transfer can be determined by mea- 
suring the response of the two systems at their ends only. 



can search numerically for optimised measurement times using Eq. (|2.44p and the 
condition from Eq. (|2.43l) . 

One weakness of the scheme described here is that the times at which Bob measures 
have to be very precise, because otherwise the measurements will not be unbiased 
with respect to a and (5. This demand can be relaxed by measuring at times where 
not only the probability amplitudes are similar, but also their slope (see Fig. 12. Zfl . 
The computation of these optimal timings for a given system may be complicated, but 
they only need to be done once. 



2.10 Numerical Examples 

In this section, we show some numerical examples for two chains with Heisenberg 
couplings J which are fluctuating. The Hamiltonians of the chains 2 = 1,2 are given 
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by 

N-l 

H® = Yl J (! + ( X i t)X nl + y^ii + Z^Z%) , (2.49) 

n=l ^ 

where 5$ are uniformly distributed random numbers from the interval [—A, A] . We 

(i) 

have considered two different cases: in the first case, the 5 n are completely uncorre- 
cted (i.e. independent for both chains and all sites along the chain). In the second 

(i) 

case, we have taken into account a spacial correlation of the signs of the 5n along each 
of the chains, while still keeping the two chains uncorrelated. For both cases, we find 
that arbitrarily perfect transfer remains possible except for some very rare realisations 
of the 5 { n ] . 

Because measurements must only be taken at times which fulfil the condition fl2.43j) . 
and these times usually do not coincide with the optimal probability of finding an 
excitation at the ends of the chains, it is clear that the probability of failure at each 
measurement will in average be higher than for chains without fluctuations. Therefore, 
more measurements have to be performed in order to achieve the same probability of 
success. The price for noisy couplings is thus a longer transmission time and a higher 
number of gating operations at the receiving end of the chains. Some averaged values 
are given in Table |2~T1 for the Heisenberg chain with uncorrelated coupling fluctuations. 





A = 


A = 0.01 


A = 0.03 


A = 0.05 


A = 0.1 




377 


524 ± 27 


694 ± 32 


775 ± 40 


1106 ±248 


M 


28 


43 ±3 


58 ±3 


65 ±4 


110 ±25 



Table 2.1: The total time t and the number of measurements M needed to achieve a 
probability of success of 99% for different fluctuation strengths A (uncorrelated case). 
Given is the statistical mean and the standard deviation. The length of the chain is 
N = 20 and the number of random samples is 10. For strong fluctuations A = 0.1, 
we also found particular samples where the success probability could not be achieved 
within the time range searched by the algorithm. 

(i) 

For the case where the signs of the 5 n are correlated, we have used the same model 
as in [99], introducing the parameter c such that 

> with propability c, (2.50) 

and 

S n )s n-l < with propability 1 - c. (2.51) 



52 



2 Dual Rail encoding 



For c = 1 (c = 0) this corresponds to the case where the signs of the couplings 
are completely correlated (anti-correlated). For c = 0.5 one recovers the case of 
uncorrected couplings. We can see from the numerical results in Table 12.21 that 
arbitrarily perfect transfer is possible for the whole range of c. 





c = 


c = 0.1 


c = 0.3 


c = 0.7 


c = 0.9 


c= 1 




666 ± 20 


725 ± 32 


755 ± 41 


797 ± 35 


882 ± 83 


714 ±41 


M 


256 ±2 


62 ±3 


65 ±4 


67 ±4 


77 ±7 


60 ±4 



Table 2.2: The total time t and the number of measurements M needed to achieve a 
probability of success of 99% for different correlations c between the couplings [see Eq. 
(I2.50P and Eq. (|2.51h ]. Given is the statistical mean and the standard deviation for a 
fluctuation strength of A = 0.05. The length of the chain is N = 20 and the number 
of random samples is 20. 

For A = 0, we know from Section l2~4l that the time to transfer a state with proba- 
bility of failure P scales as 



If we want to obtain a similar formula in the presence of noise, we can perform a fit 
to the exact numerical data. For uncorrelated fluctuations of A = 0.05, this is shown 
in Fig. 12.91 The best fit is given by 



We conclude that weak fluctuations (say up to 5%) in the coupling strengths do not 
deteriorate the performance of our scheme much for the chain lengths considered. Both 
the transmission time and the number of measurements raise, but still in a reasonable 
way [cf. Table |2~H and Fig. 12.9] . For larger fluctuations, the scheme is still applicable in 
principle, but the amount of junk (i.e. chains not capable of arbitrary perfect transfer) 
may get too large. 

Note that we have considered the case where the fluctuations 5 l n are constant in 
time. This is a reasonable assumption if the dynamic fluctuations (e.g. those arising 
from thermal noise) can be neglected with respect to the constant fluctuations (e.g. 
those arising from manufacturing errors). If the fluctuations were varying with time, 
the tomography measurements in Sec. 12.91 would involve a time-average, and Bob 
would not measure exactly at the correct times. The transferred state (|2.33|) would 
then be affected by both phase and amplitude noise. 



t(P) = 0.33J _1 iV i e) |lnP 



(2.52) 



t(P) = 0.2J~\N' 1 - 9 |lnP 



(2.53) 
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Figure 2.9: Time t needed to transfer a state with a given joint probability of failure 
P across a chain of length iV with uncorrelated fluctuations of A = 0.05. The points 
denote numerical data averaged over 100 realisations, and the fit is given by Eq. (|2.53p . 
This figure should be compared with Fig. 12.41 where A = 0. 



2.11 Coupled chains 

Let us look at the condition for conclusive transfer in the more general scenario indi- 
cated by Fig. 12.101 Alice and Bob have a black box acting as an amplitude damping 
channel in the following way. It has two inputs and two outputs. If Alice puts in state 
in the dual rail, 

|V?) = a|01) +0|1O), (2.54) 

where a and (3 are arbitrary and unknown normalised amplitudes, then the output at 
Bob is given by 

j#W| + (l-p)|00>(00|, (2.55) 
with a normalised "success" state 

\<f>) = 4= \ a f\° l ) + 09\W) + a/| 10) + (3g\01)} • (2-56) 
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This black box describes the behaviour of an arbitrarily coupled qubit system that 
conserves the number of excitations and that is initialised in the all zero state, including 
parallel uncoupled chains, and coupled chains. 



Ai- 
A 2 - 



-Bi 

I- B 2 



Figure 2.10: Most general setting for conclusive transfer: A black box with two inputs 
and two outputs, acting as an amplitude damping channel defined by Eqs. (|2.54j) and 
(l2~55l 



From the normalisation of \6) it follows that 



p = p(a,f3) = \af + pg\ + 



Pg + af 



(2.57) 



We are interested in conclusive transfer: by measuring the observable |00)(00| the 
Bob can project the output onto either the failure state 1 00) or \(j>). This is clearly 
possible, but the question is if the output and the input \tp) are related by a 
unitary operation. 

If Bob is able to recover the full information that Alice sent, then p(a,/3) must be 
independent of a and [3 (otherwise, some information on these amplitudes could be 
obtained by the measurement already, which contradicts the non-cloning theorem [2]). 
This implies that p(l, 0) = p(0, 1), i.e. 



m 2 + 



/ 



\9\ 



\9\ 



Because 



it also implies that 



= p(l,Q)+Be{rg + gf*} 
Re {/*§ + <?/*} 



0. 



(2.58) 

(2.59) 
(2.60) 

(2.61) 
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Using the same trick for p ("^j, -^j we get that Im + = and therefore 

f*9 + 9f* = 0. (2.62) 

If we write \ip) = U\(f>) we get 

1 / f f \ 

(2.63) 

which is a unitary operator if Eq. (|2.58|) and (I2,62p hold. We thus come to the 
conclusion that conclusive transfer with the black box defined above is possible if and 
only if the probability p is independent of a and (3. It is interesting to note that a 
vertical mirror symmetry of the system does not guarantee this. A counterexample 
is sketched in Fig. 1 2 . Ill clearly the initial ("dark") state 1 01) — 1 10) does not evolve, 
whereas 1 01) + 1 10) does. Hence the probability must depend on a and f3. A trivial 
case where conclusive transfer works is given by two uncoupled chains, at times where 
|/| 2 = \g\ 2 . This was discussed in Sect. 12.81 A non-trivial example is given by the 
coupled system sketched in Fig. 12.121 This can be seen by splitting the Hamiltonian 
in a horizontal and vertical component, 

H = H V + H z . (2.64) 

By applying H V H Z and H Z H V on single-excitation states it is easily checked that they 
commute in the first excitation sector (this is not longer true in higher sectors). Since 
the probability is independent of a and (3 in the uncoupled case it must also be true 
in the coupled case (a rotation in the subspace {|01), |10)} does not harm). 




Figure 2.11: A simple counterexample for a vertically symmetric system where dual 
rail encoding is not possible. The black lines represent exchange couplings. 

A final remark - as Alice and Bob alway only deal with the states {|00), 1 10) , |01)} 
it is obvious that the encoding used in this chapter is really living on qutrits. In some 
sense it would be more natural to consider permanently coupled systems of qutrits, 
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Figure 2.12: An example for a vertically symmetric system where dual rail encoding 
is possible. The black lines represent exchange couplings of equal strength. 

such as SU(3) chains [102, 123-125]. The first level of the qutrit |0) is then used as a 
marker for "no information here", whereas the information is encoded in the states |1) 
and 1 2). One would have to ensure that there is no transition between |0) and |1), |2), 
and that the system is initialised in the all zero state. 

2.12 Conclusion 

In conclusion, we have presented a simple scheme for conclusive and arbitrarily per- 
fect quantum state transfer. To achieve this, two parallel spin chains (individually 
amplitude damping channels) have been used as one amplitude delaying channel. We 
have shown that our scheme is more robust to decoherence and imperfect timing than 
the single chain schemes. We have also shown that the scheme is applicable to dis- 
ordered and coupled chains. The scheme can be used as a way of improving any of 
the other schemes from the introduction. For instance, one may try to engineer the 
couplings to have a very high probability of success already at the first measurement, 
and use further measurements to compensate the errors of implementing the correct 
values for the couplings. We remark that the dual rail protocol is unrelated to error 
filtration [126] where parallel channels are used for filtering out environmental effects 
on flying qubits, whereas the purpose of the dual rail protocol is to ensure the ar- 
rival of the qubit. Indeed one could combine both protocols to send a qubit on say 
four rails to ensure the arrival and filter errors. Finally, we note that in some recent 
work [80] it was shown that our encoding can be used to perform quantum gates while 
the state is transferred, and that it can increase the convergence speed if one performs 
measurements at intermediate positions [110,127]. 
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3.1 Introduction 

In quantum information theory the rate R of transferred qubits per channel is an 
important efficiency parameter [70]. Therefore one question that naturally arises is 
whether or not there is any special meaning in the 1/2 value of R achieved in the 
dual rail protocol of the last chapter. We will show now that this is not the case, 
because there is a way of bringing R arbitrarily close to 1 by considering multi rail 
encodings. Furthermore, in Section 12.31 it was still left open for which Hamiltonians 
the probability of success can be made arbitrarily close to 1. Here, we give a sufficient 
and easily attainable condition for achieving this goal. 

This chapter is organised as follows: the model and the notation are introduced in 
Sec. 13.21 The efficiency and the fidelity of the protocol are discussed in Sec. 13.31 and 
in Sec. 13.41 respectively. Finally in Sec. 13.51 we prove a theorem which provides us 
with a sufficient condition for achieving efficient and perfect state transfer in quantum 
chains. 

3.2 The model 

Assume that the two communicating parties operate on M independent (i.e. non 
interacting) copies of the chain. This is quite a common attitude in quantum informa- 
tion theory [70] where successive uses of a memoryless channel are formally described 
by introducing many parallel copies of the channel (see [54] for a discussion on the 
possibility of applying this formal description to quantum chain models). Moreover for 
the case at hand the assumption of Alice and Bob dealing with "real" parallel chains 
seems reasonable also from a practical point of view [48,49]. The idea is to use these 
copies to improve the overall fidelity of the communication. As usual, we assume Alice 
and Bob to control respectively the first and last qubit of each chain (see Fig. 13. lj) . By 
preparing any superposition of her spins Alice can in principle transfer up to M logical 
qubits. However, in order to improve the communication fidelity the two parties will 
find it more convenient to redundantly encode only a small number (say Q{M) ^ M) 
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Length N 





Figure 3.1: Schematic of the system: Alice and Bob operate M chains, each containing 
N spins. The spins belonging to the same chain interact through the Hamiltonian H 
which accounts for the transmission of the signal in the system. Spins of different 
chains do not interact. Alice encodes the information in the first spins of the chains 
by applying unitary transformations to her qubits. Bob recovers the message in the 
last spins of the chains by performing joint measurements. 

of logical qubits in the M spins. By adopting these strategies Alice and Bob are ef- 
fectively sacrificing the efficiency R(M) = Q(M)/M of their communication line in 
order to increase its fidelity. This is typical of any communication scheme and it is 
analogous to what happens in quantum error correction theory, where a single logical 
qubit is stored in many physical qubits. In the last chapter we have seen that for 
M = 2 it is possible to achieve perfect state transfer of a single logical qubit with an 
efficiency equal to 1/2. Here we will generalise such result by proving that there exist 
an optimal encoding-decoding strategy which asymptotically allows to achieve perfect 
state transfer and optimal efficiency, i.e. 

lim R(M) = 1 . (3.1) 

Our strategy requires Alice to prepare superpositions of the M chains where ~ M/2 
of them have a single excitation in the first location while the remaining are in |0). 
Since in the limit M » 1 the number of qubit transmitted is log (j^y 2 ) ~ M, this 
architecture guarantees optimal efficiency (|3.ip . On the other hand, our protocol 
requires Bob to perform collective measurements on his spins to determine if all the 
~ M/2 excitations Alice is transmitting arrived at his location. We will prove that 
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by repeating these detections many times, Bob is able to recover the messages with 
asymptotically perfect fidelity. 

Before beginning the analysis let us introduce some notation. The following defi- 
nitions look more complicated than they really are; unfortunately we need them to 
carefully define the states that Alice uses for encoding the information. In order to 
distinguish the M different chains we introduce the label m = 1, • • • , M: in this for- 
malism \n) m represents the state of m-th chain with a single excitation in the n-th 
spin. In the following we will be interested in those configurations of the whole system 
where K chains have a single excitation while the remaining M — K are in |0), as in 
the case 

|l)i ® |1> 2 ■ ■ ■ ® 8> |0>jf+i • • ■ ® \0) M (3.2) 

where for instance the first K chains have an excitation in the first chain location. 
Another more general example is given in Fig. 13.21 The complete characterisation of 
these vectors is obtained by specifying i) which chains possess a single excitation and 
%%) where these excitations are located horizontally along the chains. In answering to 
the point i) we introduce the if-element subsets Se, composed by the labels of those 
chains that contain an excitation. Each of these subsets Si corresponds to a subspace 
of the Hilbert space TL(Si) with a dimension N K . The total number of such subsets 
is equal to the binomial coefficient (^) , which counts the number of possibilities in 
which K objects (excitations) can be distributed among M parties (parallel chains). 
In particular for any I = 1, • • • , ( K ) the £-th subset Si will be specified by assigning 
its K elements, i.e. S e = {mf , ■ ■ ■ , mj?} with mf G {1, • • • , M} and mf < m% 
for all j = 1, • • • , K. To characterise the location of the excitations, point ii), we will 
introduce instead the K-dimensional vectors n = (m, • • • , uk) where rij G {1, • • ■ , N}. 
We can then define 

K 

\n;£)) = (g) \ nj ) m w (g) |0> m , , (3.3) 

i =1 3 m'eSe 

where Si is the complementary of Si to the whole set of chains. 

The state (|3.3|) represents a configuration where the j-th chain of the subset Si is in 
\rij) while the chains that do not belong to Si are in |0) (see Fig. 13.21 for an explicit 
example). The kets \n;£)) are a natural generalisation of the states \n)i 1 0)2 and 
|0)i<g>|n}2 which were used for the dual rail encoding. They are useful for our purposes 
because they are mutually orthogonal, i.e. 

((n;l\n , ;£')) = 6 w S m> , (3.4) 
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Length 6 



a 

in 




isXXXXX) 




A 



B 



Figure 3.2: Example of our notation for M = 5 chains of length N = 6 with if = 2 
excitations. The state above, given by |0)i &> 1 3)2 ® 1 0)3 (g) 1 1)4 (8> 1 0)5 , has excitations 
in the chains mi = 2 and 771,2 = 4 at the horizontal position n\ = 3 and 712 = 1. It 
is in the Hilbert space T-1{Sq) corresponding to the subset Sq = {2,4} (assuming that 
the sets Se are ordered in a canonical way, i.e. Si = {1, 2}, S2 = {1, 3} and so on) and 
will be written as |(3, 1);6)). There are ( 2 ) = 10 different sets Se and the number of 
qubits one can transfer using these states is log 2 10 as 3. The efficiency is thus given 
by R « 3/5 which is already bigger than in the dual rail scheme. 



and their time evolution under the Hamiltonian does not depend on I. Among the 
vectors (13. 3p those where all the K excitations are located at the beginning of the 
chains play an important role in our analysis. Here n = 1 = (1, • • • ,1) and we can 
write 

= |l) m (g) |0} m , . (3.5) 

™G5f m'eSi 

According to Eq. (|3.4|) . for I = 1, • • • , ( K ) these states form orthonormal set of (^) 
elements. Analogously by choosing n = N = (N, ■ ■ ■ , N) we obtain the orthonormal 
set of ( K ) vectors 

\N;£}} = (g) \N) m (g) |0) m ,, (3.6) 

m£S e m'eS e 

where all the K excitations are located at the end of the chains. 



61 



3 Multi Rail encoding 



3.3 Efficient encoding 

If all the M chains of the system are originally in |0), the vectors (13.51) can be prepared 
by Alice by locally operating on her spins. Moreover since these vectors span a ( K ) 
dimensional subspace, Alice can encode in the chain Q(M, K) = log 2 ( K ) qubits of 
logical information by preparing the superpositions, 

\$))=J2Ae\i;l)), (3.7) 

i 

with An complex coefficients. The efficiency of such encoding is hence R(M, K) = 
l ° S2 M K ^ which maximised with respect to K gives, 

1 ( !og 2 (m/o) for M even 

R(M) = —{ 2 ) M/ Z . (3.8) 
1 log 2 ( ( m M i)/ 2 ) for M odd . 

The Stirling approximation can then be used to prove that this encoding is asymptot- 
ically efficient (|3.ip in the limit of large M, e.g. 



/ M \ M M 
\M/2J ~ g2 (M/2V 



fog 2 (, /f/ J « l og 27T7 — M = M. (3.9) 



Note that already for M = 5 the encoding is more efficient (cf. Fig. I3.2p than in the 
dual rail encoding. In the remaining of the chapter we show that the encoding (|3.7I) 
provides perfect state transfer by allowing Bob to perform joint measurements at his 
end of the chains. 



3.4 Perfect transfer 

Since the M chains do not interact with each other and possess the same free Hamil- 
tonian H, the unitary evolution of the whole system is described by U{t) = ® m ii m (t), 
with u m (t) being the operator acting on the m-th chain. The time evolved of the input 
1 1; £)) of Eq. J3J2) is thus equal to 

U(t)\l;£}} = ^F[n,l;t] \n;£)) , (3.10) 

ft 

where the sum is performed for all rij = 1, • • ■ , N and 

F[n,n';t] = f nu< (t) ■ ■ ■ f nx , n ' K (t) , (3.11) 
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is a quantity which does not depend on I. In Eq. (13 . 10|) the term n = N corresponds 
to having all the K excitations in the last locations of the chains. We can thus write 



u(t)\T-e)) = 7i(*)l^» + \A - l7i(*)l 2 , (3.12) 

where 

7 i(t) = ((N;l\U(t)\l;£))=F[N,l;t] (3.13) 

is the probability amplitude that all the K excitation of arrive at the end of the 
chains, and 

m);Q= ^Fi[n,l;i]|n;£)), (3.14) 

with 

F 1 [n,j;t] = ^ [n,T;t] , (3.15) 
Vl-I7i(i)| 2 

is a superposition of terms where the number of excitations arrived to the end of the 
communication line is strictly less then K. It is worth noticing that Eq. (13. 4ft yields 
the following relations, 

«#;4£(t);O=0, ((memY,e'))=5ti' , (3-16) 

which shows that 1S an orthonormal set of vectors which spans a subspace 

orthogonal to the states \N;£)). The time evolution of the input state (|3.7p follows by 
linearity from Eq. (|3.12l) . i.e. 



with 



|*(t)» = 7l(t) I*)) + Vl-lTiWI 2 |*(f)» , (3.17) 
|*» ee (3.18) 



The vectors |\&)) and \^(t))) are unitary transformations of the input message (|3.7|) 
where the orthonormal set {|1;^))} has been rotated into {\N;£)}} and {\£(t);£}}} 
respectively. Moreover |^)) is the configuration we need to have for perfect state 
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transfer at the end of the chain. In fact it is obtained from the input message (|3.7I) 
by replacing the components |1) (excitation in the first spin) with \N) (excitation in 
the last spin). From Eq. (|3.16|) we know that |^)) and (^(i))) are orthogonal. This 
property helps Bob to recover the message \*S>)) from |<3?(i))): he only needs to perform 
a collective measurement on the M spins he is controlling to establish if there are K 
or less excitations in those locations. The above is clearly a projective measurement 
that can be performed without destroying the quantum coherence associated with the 
coefficients A(_. Formally this can described by introducing the observable 

0= l-£ |#;*)>((#; *| ■ ( 3 ' 19 ) 

t 

A single measurement of 6 on |<3?(ii))) yields the outcome with probability p\ = 
|7i(ti)| 2 , and the outcome +1 with probability 1 — p\. In the first case the system 
will be projected in |^)) and Bob will get the message. In the second case instead the 
state of the system will become |\&(£i))). Already at this stage the two communicating 
parties have a success probability equal to p\. Moreover, as in the dual rail protocol, 
the channels have been transformed into a quantum erasure channel [113] where the 
receiver knows if the transfer was successful. Just like the dual rail encoding, this 
encoding can be used as a simple entanglement purification method in quantum chain 
transfer (see end of Section l2~2j) . The rate of entanglement that can be distilled is 
given by 

R(M) F[N, 1; t] 2 = R{M)p{t) L M / 2 J , (3.20) 

where we used Eq. (|3.1ip and p(t) = |/jv,i(i)| 2 . As we can see, increasing M on 
one hand increases R(M), but on the other hand decreases the factor p(t)L M / 2 J. Its 
maximum with respect to M gives us a lower bound of the entanglement of distillation 
for a single spin chain, as shown in Fig. 11.111 We can also see that it becomes worth 
encoding on more than three chains for conclusive transfer only when p(t) > 0.8. 

Consider now what happens when Bob fails to get the right answer from the mea- 
surement. The state on which the chains is projected is explicitly given by 

= J2F 1 [n,l;ti]J2Mn;t)) ■ (3.21) 

Let us now consider the evolution of this state for another time interval t2- By repeat- 
ing the same analysis given above we obtain an expression similar to (I3.17p . i.e. 



I*(t2,ti)» = 72 |*)) + Vl-|72| 2 |*(t 2 ,tx)» , (3.22) 
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where now the probability amplitude of getting all excitation in the N-th locations is 
described by 

72 = Y / F[N,n;t 2 ]F 1 [n,l;ti]. (3.23) 

In this case \^(t)}} is replaced by 

|*(*2,ti)» = J2A e \Z(t 2 ,ti);£)) , (3.24) 

e 

with 

m2,h);£)) = ]T F 2 [n,l;t 2 M]\n-A)), (3.25) 

and F 2 defined as in Eq. (|3.27p (see below). In other words, the state \&(t 2 , h))} can 
be obtained from Eq. (|3.17|) by replacing 71 and F\ with 72 and F 2 . Bob can hence 
try to use the same strategy he used at time t±: i.e. he will check whether or not 
his M qubits contain K excitations. With (conditional) probability p 2 = \~/ 2 \ 2 he will 
get a positive answer and his quantum register will be projected in the state \^f}} of 
Eq. (I3.18p . Otherwise he will let the system evolve for another time interval £3 and 
repeat the protocol. By reiterating the above analysis it is possible to give a recursive 
expression for the conditional probability of success p q = \^ q \ 2 after q — 1 successive 
unsuccessful steps. The quantity 7^ is the analogue of 72 and 71 of Eqs. (|3.13|) and 
(I3~22l . It is given by 

lq = F[N,n;t q ] F q ^[H, 1, t q -i, ■ ■ ■ ,h] , (3.26) 

where 

F,_i[n,l; M] (3.27) 

and Fi[n, l,t] is given by Eq. (I3.15p . In these equations t q , ■ ■ ■ ,ti are the time-intervals 
that occurred between the various protocol steps. Analogously the conditional proba- 
bility of failure at the step q is equal to 1 — p q . The probability of having j — 1 failures 
and a success at the step j-th can thus be expressed as 

<j) = Pi(l-Pj-l)(l-Pj-2) "•(!- Pi), (3.28) 
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while the total probability of success after q steps is obtained by the sum of n(j) for 
all j = !,■■■ ,<?, i-e- 



Since pj ^ 0, Eq. (|3.29p is a monotonic function of q. As a matter of fact in the next 
section we prove that under a very general hypothesis on the system Hamiltonian, the 
probability of success P q converges to 1 in the limit of q — > oo. This means that by 
repeating many times the collective measure described by B Bob is guaranteed to get, 
sooner or later, the answer and hence the message Alice sent to him. In other words 
our protocol allows perfect state transfer in the limit of repetitive collective measures. 
Notice that the above analysis applies for all classes of subsets Sg. The only difference 
between different choices of K is in the velocity of the convergence of P q — > 1. In any 
case, by choosing K ~ M/2 Alice and Bob can achieve perfect fidelity and optimal 
efficiency. 

3.5 Convergence theorem 



Theorem 3.1 (Arbitrarly perfect transfer) If there is no eigenvector \e m ) of 
the quantum chain Hamiltonian H which is orthogonal to \N) , then there is a 
choice of the times intervals t q ,t q -i, ■ ■ ■ ,ti such that the fidelity converges to 1 as 
q — ► oo. 



Before proving this Theorem, let us give an intuitive reasoning for the convergence. 
The unitary evolution can be thought of of a rotation in some abstract space, while 
the measurement corresponds to a projection. The dynamics of the system is then 
represented by alternating rotations and projections. In general this will decrease the 
norm of each vector to null, unless the rotation axis is the same as the projection axis. 

Proof The state of the system at a time interval of t q after the (q — l)-th failure can 
be expressed in compact form as follows 



with U(t) the unitary time evolution generated by the system Hamiltonian, and with 




(3.29) 



u(t q )eu(t q ^ 1 )Q---u(t 1 )@m 

y/(l- Pq -l)---(l-pi) 



(3.30) 
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the projection denned in Eq. (|3.19j) . One can verify for instance that for q = 2, the 
above equation coincides with Eq. (|3.22p . [For q = 1 this is just (|3.17|) evaluated at 
time t±]. By definition the conditional probability of success at step g-th is equal to 

p, = |((*|*(t„... , *i)»| 2 - (3.31) 

Therefore, Eq. (EHSl) yields 

n{q) = |((^|C/(t 9 )e[/(Vi)©---^(ii)©l^))| 2 (3-32) 
= \((N;£\U(t g )eU(t q ^)G---U(tt)e\l;£))f , 

where the second identity stems from the fact that, according to Eq. (|3.4p . U(t)@ 
preserves the orthogonality relation among states \n;£)) with distinct values of I. In 
analogy to the cases of Eqs. (13. lip and (I3.13p . the second identity of (|3.32j) establishes 
that 7r(g) can be computed by considering the transfer of the input for arbitrary 

I. The expression (|3.32|) can be further simplified by noticing that for a given I 
the chains of the subset Se contribute with a unitary factor to ir(q) and can be thus 
neglected (according to (13. 5p they are prepared in |0) and do not evolve under U(t)&). 
Identify \l))e and \N))e with the components of \1;£)) and \N;£)) relative to the chains 
belonging to the subset Se. In this notation we can rewrite Eq. (I3.32p as 

7r(q) = u(N\u e (t q )e e ■■■u e (t 1 )e e \i)} e \ 2 , (3.33) 

where @e = 1 — |JV))^((JV| and Ue(t) is the unitary operator <Xi m g5 £ n m (t) which de- 
scribes the time evolution of the chains of Sc. To prove that there exist suitable 
choices of t£ such that the series (|3.29|) converges to 1 it is sufficient to consider the 
case ti = t > for all j = 1, • • • ,q: this is equivalent to selecting decoding proto- 
cols with constant measuring intervals. By introducing the operator Te = Ue(t)@e, 
Eq. (|3.33|) becomes thus 

7r(q) = U(N\ (T e y\l)) e \ 2 (3.34) 
= e ((l\(T}y\NU(N\ (T e y\l)) e = w(q) - w(q + 1) , 

where 

w(j) ^((l\(T}y (T e y\i)) e = ||(T^|r»,|| 2 , (3.35) 
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is the norm of the vector (T e y\T))t. Substituting Eq. (ET33jl in Eq. ff3T29l) yields 

P q = Y,[w(j)-w(j + l)] = l-w(q+l) (3.36) 
i=i 

where the property w(l) = ^((1|©^| 1))^ = 1 was employed. Proving the thesis is hence 
equivalent to prove that for q — > oo the succession w(q) nullifies. This last relation can 
be studied using properties of power bounded matrices [128]. In fact, by introducing 
the norm of the operator (Ti) q we have, 

w(q) = [|CZi)«|i»if < \myf < c (l±fM] 2q (3. 37) 



where c is a positive constant which does not depend on q (if S is the similarity 
transformation that puts Tg into the Jordan canonical form, i.e. J = S~ 1 TgS, then 
c is given explicitly by c = \\S\\ US' -1 !!) and where p(Tt) is the spectral radius of Tg, 
i.e. the eigenvalue of Ti with maximum absolute value (N.B. even when Tg is not 
diagonalisable this is a well defined quantity). Equation (|3.37|) shows that p(Tt) < 1 
is a sufficient condition for w(q) — > 0. In our case we note that, given any normalised 
eigenvector |A))^ of Tg with eigenvalue A we have 

|A| = ||2}|A»*|| = ||0*|A»,|| < 1 , (3.38) 

where the inequality follows from the fact that is a projector. Notice that in 
Eq. (|3.38p the identity holds only if |A)) is also an eigenvector of 0^ with eigenvalue +1, 
i.e. only if |A))^ is orthogonal to \N))g. By definition |A))^ is eigenvector Tg = Ui(t)&£: 
therefore the only possibility to have the equality in Eq. (|3.38 



lianj. 



is that i) is 

an eigenvector of Ue(t) (i.e. an eigenvector of the Hamiltoniarlll Hg ^ of the chain 
subset St) and ii) it is orthogonal to |iV))^. By negating the above statement we 
get a sufficient condition for the thesis. Namely, if all the eigenvectors \E))i of H\ ^ 
are not orthogonal to \N))t than the absolute values of the eigenvalues A of Tg are 
strictly smaller than 1 which implies piTg) < 1 and hence the thesis. Since the Sg 
channels are identical and do not interact, the eigenvectors \E))t = <S> m( zs e \ e m)m are 
tensor product of eigenvectors \e m ) of the single chain Hamiltonians H . Therefore the 



1 Notice that strictly speaking the eigenvectors of the Hamiltonian are not the same as those of the 
time evolution operators. The latter still can have evolution times at which additional degeneracy 
can increase the set of eigenstates. A trivial example is given for t = where all states become 
eigenstates. But it is always possible to find times t at which the eigenstates of U(t) coincide with 
those of H . 
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sufficient condition becomes 

e((E\N)) e = J[ m (N\e m ) m + , (3.39) 

m£Si 

which can be satisfied only if (N\e m ) ^ for all eigenvectors \e m ) of the single chain 
Hamiltonian H. ■ 

Remark 3.1 While we have proven here that for equal time intervals the probability 
of success is converging to unity, in practice one may use optimal measuring time 
intervals U for a faster transfer (see also Section f2T4|) . We also point out that timing 
errors may delay the transfer, but will not decrease its fidelity. 



3.6 Quantum chains with nearest-neighbour interactions 

It is worth noticing that Eq. (13.39P is a very weak condition, because eigenstates 
of Hamiltonians are typically entangled. For instance, it holds for open chains with 
nearest neighbour-interactions: 



Theorem 3.2 (Multi rail protocol) Let H be the Hamiltonian of an open 
nearest-neighbour quantum chain that conserves the number of excitations. If there 
is a time t such that /i,jv(£) 7^ (i.e. the Hamiltonian is capable of transport be- 
tween Alice and Bob) then the state transfer can be made arbitrarily perfect by 
using the multi rail protocol. 



Proof We show by contradiction that the criterion of Theorem 13.11 is fulfilled. As- 
sume there exists a normalised eigenvector |e) of the single chain Hamiltonian H such 
that 

(N\e) = 0. (3.40) 
Because |e) is an eigenstate, we can conclude that also 

(e \H\ N) = 0. (3.41) 

If we act with the Hamiltonian on the ket in Eq. (I3.4ip we may get some term propor- 
tional to (e|iV) (corresponding to an Ising-like interaction) and some part proportional 
to (e\N — 1) (corresponding to a hopping term; if this term did not exist, then clearly 
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fi,N(t) = for all times). We can thus conclude that 

(e|JV-l) = 0. (3.42) 

Note that for a closed chain, e.g. a ring, this need not be the case, because then also 
a term proportional to {e\N + 1) = (e|l) would occur. If we insert the Hamiltonian 
into Eq. (j3.42|) again, we can use the same reasoning to see that 

(e\N — 2) = • • • = (e|l) = (3.43) 

and hence |e) = 0, which is a contradiction to |e) being normalised. ■ 



3.7 Comparison with Dual Rail 

As we have seen above, the Multi Rail protocol allows us in principle to reach in 
principle a rate arbitrarily close to one. However for a fair comparison with the Dual 
Rail protocol, we should also take into account the time-scale of the transfer. For the 
conclusive transfer of entanglement, we have seen in Section I3~4l that only for chains 
which have a success probability higher than p(t) = 0.8 it is worth encoding on more 
than three rails. The reason is that if the probability of success for a single excitation 
is p, then the probability of success for |_-^/2j excitations on on M parallel chains is 
lowered to pL M / 2 J . The protocol for three rails is always more efficient than on two, as 
still only one excitation is being used, but three complex amplitudes can be transferred 
per usage. 

For arbitrarily perfect transfer, the situation is slightly more complicated as the 
optimal choice of M also depends on the joint probability of failure that one plans to 
achieve. Let us assume that at each step of the protocol, the success probability on a 
single chain is p. Then the number of steps to achieve a given probability of failure P 
using M chains is given by 

£(P,M) = max{ in(i ^ M/2j) ,l}. (3.44) 

If we assume that the total time-scale of the transfer is proportional to the number of 
steps, then the number of qubits that can be transferred per time interval is given by 

v(P,M) oc R(M)/l{P,M). (3.45) 

Optimising this rate with respect to M we find three different regimes of the joint 
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probability of failure (see Fig. 13. 3h , If one is happy with a large P, then the Multi Rail 
protocol becomes superior to the Dual Rail for medium p. For intermediate P, the 
threshold is comparable to the threshold of p = 0.8 for conclusive transfer of entangle- 
ment. Finally for very low P the Multi Rail only becomes useful for p very close to one. 
In all three cases the threshold is higher than the p(t) that can usually achieved with 
unmodulated Heisenberg chains. We can thus conclude that the Multi Rail protocol 
only becomes useful for chains which already have a very good performance. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Success probability on a single chain 



Figure 3.3: Optimal rates (maximisation of Eq. (I3.45l with respect to M) for the Multi 
Rail protocol. Shown are three curves corresponding to different values of the joint 
probability of failure P one plans to achieve. 



3.8 Conclusion 

We thus conclude that any nearest-neighbour Hamiltonian that can transfer quantum 
information with nonzero fidelity (including the Heisenberg chains analysed above) is 
capable of efficient and perfect transfer when used in the context of parallel chains. 
Hamiltonians with non-nearest neighbour interactions [81,89] can also be used as long 
as the criterion of Theorem 13.11 is fulfilled. 
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4.1 Introduction 

We have seen above that by applying measurements at the end of parallel chains, the 
state of the chain is converging to the ground state, and the quantum information is 
transferred to the receiver. Indeed, repetitive application of the same transformation 
is the key ingredient of many controls techniques. Beside quantum state transfer, they 
have been exploited to inhibit the decoherence of a system by frequently perturbing 
its dynamical evolution [129-133] (Bang-Bang control) or to improve the fidelity of 
quantum gates [134] by means of frequent measurements (quantum Zeno-effect [135]). 
Recently analogous strategies have also been proposed in the context of state prepa- 
ration [136-142]. In Refs. [138,139] for instance, a homogenisation protocol was pre- 
sented which allows one to transform any input state of a qubit into a some pre-fixed 
target state by repetitively coupling it with an external bath. A similar thermalisa- 
tion protocol was discussed in Ref. [140] to study the efficiency of simulating classical 
equilibration processes on a quantum computer. In Refs. [141,142] repetitive interac- 
tions with an externally monitored environment were exploited instead to implement 
purification schemes which would allow one to extract pure state components from 
arbitrary mixed inputs. 




ergodic mixing 



Figure 4.1: Schematic examples of the orbits of a ergodic and a mixing map. 



72 



4 Ergodicity and mixing 



The common trait of the proposals [136-142] and the dual and multi rail protocols 
is the requirement that repeated applications of a properly chosen quantum operation 
r converges to a fixed density matrix independently from the input state x of the 
system, i.e. 

— ► x* , (4.1) 



with "o" representing the composition of maps. Following the notation of Refs. [143, 
144] we call Eq. (|4.ip the mixing property of r. It is related with another important 
property of maps, namely ergodicity (see Fig. 14. ip . The latter requires the existence of 
a unique input state xq which is left invariant under a single application of the map0, 
i.e. 

t(x) = x <J=^ x = Xq . (4.2) 

Ergodicity and the mixing property are of high interest not only in the context of the 
above quantum information schemes. They also occur on a more fundamental level 
in statistical mechanics [147] and open quantum systems [121,148], where one would 
like to study irreversibility and relaxation to thermal equilibrium. 

In the case of quantum transformations one can show that mixing maps with con- 
vergence point x* are also ergodic with fixed point xq = x*. The opposite implication 
however is not generally true since there are examples of ergodic quantum maps which 
are not mixing (see the following). Sufficient conditions for mixing have been discussed 
both in the specific case of quantum channel [140, 143, 146] and in the more abstract 
case of maps operating on topological spaces [147]. In particular the Lyapunov direct 
method [147] allows one to prove that an ergodic map r is mixing if there exists a con- 
tinuous functional S which, for all points but the fixed one, is strictly increasing under 
t. Here we strengthen this criterion by weakening the requirement on S: our gener- 
alised Lyapunov functions are requested only to have limiting values S(T n (x))\ n ^ oc 
which differ from S(x) for all x ^ xq. It turns out that the existence of such S is 
not just a sufficient condition but also a necessary condition for mixing. Exploiting 
this fact one can easily generalise a previous result on strictly contractive maps [143] 



1 Definition l14.2[| may sound unusual for readers who are familiar with a definition of ergodicity from 
statistical mechanics, where a map is ergodic if its invariant sets have measure or 1. The notion 
of ergodicity used here is completely different, and was introduced in [143, 145, 146]. The set X 
one should have in mind here is not a measurable space, but the compact convex set of quantum 
states. A perhaps more intuitive definition of ergodicity based on the time average of observables 
is given by Lemma l4.5p . 
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by showing that maps which are asymptotic deformations (see Definition I4,14|) are 
mixing. This has, unlike contractivity, the advantage of being a property independent 
of the choice of metric (see however [144] for methods of finding "tight" norms). In 
some cases, the generalised Lyapunov method permits also to derive an optimal mix- 
ing condition for quantum channels based on the quantum relative entropy Finally 
a slightly modified version of our approach which employs multi-central Lyapunov 
functions yields a characterisation of (not necessarily mixing) maps which in the limit 
of infinitely many applications move all points toward a proper subset (rather than a 
single point) of the input space. 

The introduction of a generalised Lyapunov method seems to be sound not only 
from a mathematical point of view, but also from a physical point of view. In effect, 
it often happens that the informations available on the dynamics of a system are only 
those related on its asymptotic behaviour (e.g. its thermalisation process), its finite 
time evolution being instead difficult to characterise. Since our method is explicitly 
constructed to exploit asymptotic features of the mapping, it provides a more effective 
way to probe the mixing property of the process. 

Presenting our results we will not restrict ourself to the case of quantum operations. 
Instead, following [147] we will derive them in the more general context of continuous 
maps operating on topological spaces [149]. This approach makes our results stronger 
by allowing us to invoke only those hypotheses which, to our knowledge, are strictly 
necessary for the derivation. It is important to stress however that, as a particular 
instance, all the Theorems and Lemmas presented in this chapter hold for any linear, 
completely positive, trace preserving map (i.e. quantum channels) operating on a 
compact subset of normed vectors (i.e. the space of the density matrices of a finite 
dimensional quantum system). Therefore readers who are not familiar with topological 
spaces can simply interpret our derivations as if they were just obtained for quantum 
channels acting on a finite dimensional quantum system. 

This chapter is organised as follows. In Sec. 14.31 the generalised Lyapunov method 
along with some minor results is presented in the context of topological spaces. Then 
quantum channels are analysed in Sec. 14.41 providing a comprehensive summary of the 
necessary and sufficient conditions for the mixing property of these maps. Conclusions 
and remarks form the end of the chapter in Sec. 14.51 

4.2 Topological background 

Let us first introduce some basic topological background required for this chapter. A 
more detailed introduction is given in [149]. Topological spaces are a very elegant way 
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of defining compactness, convergence and continuity without requiring more than the 
following structure: 

Definition 4.1 A topological space is a pair (X, O) of a set X and a set O of subsets 
of X (called open sets) such that 

1. X and are open 

2. Arbitrary unions of open sets are open 

3. Intersections of two open sets are open 

Example 4-1 If X is an arbitrary set, and O = {Af, 0}, then (X,0) is a topological 
space. O is called the trivial topology. 

Definition 4.2 A topological space X is compact if any open cover (i.e. a set of open 
sets such that X is contained in their union) contains a finite sub-cover. 

Definition 4.3 A sequence x n 6 X is convergent with limit x* if each open neigh- 
bourhood 0(x*) (i.e. a set such that x* 6 0(x*) G contains all but finitely many 
points of the sequence. 

Definition 4.4 A map on a topological space is continous if the preimage of any open 
set is open. 

This is already all we require to make useful statements about ergodicity and mixing. 
However, there are some subtleties which we need to take care of: 

Definition 4.5 A topological space is sequentially compact if every sequence has a 
convergent subsequence. 

Sequentially compactness is in general not related to compactness! Another subtlety is 
that with the above definition, a sequence can converge to many different points. For 
example, in the trivial topology, any sequence converges to any point. This motivates 

Definition 4.6 A topological space is Hausdorff if any two distinct points can by 
separated by open neighbourhoods. 

A limit of a sequence in a Hausdorff space is unique. All these problems disappear in 
metrical spaces: 
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Definition 4.7 A metric space is a pair (X, d) of a set X and a function d : X x X — ► M 
such that 

1. c£(x, y) > and y) = <^ x = y 

2. d(x,y) = 

3. d(a;, 2) < d(x, y) + d(y, z) 

A metric space becomes a topological space with the canonical topology 

Definition 4.8 A subset O of a metric space X is open if Vx G O there is an e > 
such that {y £ X\d(x,y) <e}cO. 

In a metric space with the canonical topology, compactness and sequentially compact- 
ness become equivalent. Furthermore, it is automatically Hausdorff (see Fig. 14. 2|) . 



Figure 4.2: Relations between topological spaces [149]. The space of density matrices 
on which quantum channels are defined, is a compact and convex subset of a normed 
vectors space (the space of linear operators of the system) which, in the above graphical 
representation fits within the set of compact metric spaces. 

4.3 Generalised Lyapunov Theorem 
4.3.1 Topological spaces 

In this section we introduce the notation and derive our main result (the Generalised 
Lyapunov Theorem). 



Topological spaces 
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Definition 4.9 Let X be a topological space and let r : X —* X be a map. The 
sequence x n = T n (x), where r" is a short-hand notation for the n— fold composition 
of r, is called the orbit of x. An element x* S X is called a /jxee? point of r if and only 
if 

r(x*) = x* . (4.3) 

r is called ergodic if and only if it has exactly one fixed point, r is called mixing if 
and only if there exists a convergence point x* £ X such that any orbit converges to 
it, i.e. 

lim x n = x* \/x £ X . (4.4) 

n— >oo 

A direct connection between ergodicity and mixing can be established as follows. 

Lemma 4.1 Let t : X ^ X be a continuous mixing map on a topological Hausdorff 
space X. Then r is ergodic. 

Proof Let x* be the convergence point of r and let x £ X arbitrary. Since r is 
continuous we can perform the limit in the argument of r, i.e. 

r(x«) = r ( lim r n (x)) = lim r n+1 (x) = x*, (4.5) 

\n— >oo / n— »oo 

which shows that fixed point of r. To prove that it is unique assume by 

contradiction that r possesses a second fixed point y* ^ x*. Then linin^oo r n (y* ) = 
^ x*, so r could not be mixing (since the limit is unique in a Hausdorff space - see 
Fig. 14. 2|) . Hence r is ergodic. ■ 

Remark 4-1 The converse is not true in general, i.e. not every ergodic map is mixing 
(not even in Hausdorff topological spaces). A simple counterexample is given by 
t : [—1,1] — > [—1,1] with r(x) = —x and the usual topology of M, which is ergodic with 
fixed point 0, but not mixing since for x / 0, r n (x) = (— l) n x is alternating between 
two points. A similar counterexample will be discussed in the quantum channel section 
(see Example 14.21) . 

A well known criterion for mixing is the existence of a Lyapunov function [147]. 

Definition 4.10 Let r : X — > X be a map on a topological space X. A continuous 
map S : X — > R is called a (strict) Lyapunov function for r around x* G X if and only 
if 

S(t(x))>S{x) Vx^x*. (4.6) 
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Remark 4-2 At this point is is neither assumed that fixed point, nor that r is 

ergodic. Both follows from the theorem below. 

Theorem 4.1 (Lyapunov function) Let r : X — > X be a continuous map on a 
sequentially compact topological space X . Let S : X — > R be a Lyapunov function for 
t around x*. Then r is mixing with the fixed point x*. 

The proof of this theorem is given in [147]. We will not reproduce it here, because 
we will provide a general theorem that includes this as a special case. In fact, we will 
show that the requirement of the strict monotonicity can be much weakened, which 
motivates the following definition. 

Definition 4.11 Let r : X — > X be a map on a topological space X. A continuous 
map 5 : X — > R is called a generalised Lyapunov function for r around x* £ X if and 
only if the sequence S (r n (x)) is point-wise convergent^] for any x £ X and S fulfils 

SJx) = lim S(t h (x)) + S(x) Vx + x*. (4.7) 

In general it may be difficult to prove the point-wise convergence. However if S is 
monotonic under the action of r and the space is compact, the situation becomes 
considerably simpler. This is summarised in the following Lemma. 

Lemma 4.2 Let r : X — > X be map on a compact topological space. A continuous 
map S : X — ► R which fulfils 

S(t(x))^S(x) VxeX, (4.8) 

and 

SJx) = lim S(r n (x)) > S(x) Vx / x*. (4.9) 

for some fixed x* G X is a generalised Lyapunov function for r around x* . 

Proof It only remains to show the (point-wise) convergence of S^r^x)). Since S is 
a continuous function on a compact space, it is bounded. By Eq. (|4,8I) the sequence 
is monotonic. Any bounded monotonic sequence converges. ■ 

Corollary 4.1 Let r : X — > X be a map on a compact topological space. A continuous 
map S : X — > R which fulfils 

S(t(x))^S{x) VxGA", (4.10) 

2 Point-wise convergence in this context means that for any fixed x the sequence S n = S(r n (x)) is 
convergent. 
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and 

S(t n (x)) > S(x) Vx^x*, (4.11) 

for some fixed N € N anJ /or some x* 6 A? is a generalised Lyapunov function for r 
around x*. 

Remark 4-3 This implies that a strict Lyapunov function is a generalised Lyapunov 
function (with N = 1). 

We can now state the main result of this section: 



Theorem 4.2 (Generalized Lyapunov function) Let r : X — > X be a con- 
tinuous map on a sequentially compact topological space X. Let S : X — > R be a 
generalised Lyapunov function for r around x*. Then r is mixing with fixed point 



Proof Consider the orbit x n = T n (x) of a given x £ X. Because X is sequen- 
tially compact, the sequence x n has a convergent subsequence (see Fig. 14 .2D . i.e. 
linifc^oo x nk = x. Let us assume that x ^ x* and show that this leads to a con- 
tradiction. By Eq. (|4.7|l we know that there exists a finite N £ N such that 



S{r N (x)) jkS(x). (4.12) 
Since t n is continuous we can perform the limit in the argument, i.e. 

lim r N (x nk )=r N (x). (4.13) 

fe— >oo 

Likewise, by continuity of S we have 

lim S(x nk ) = S(x), (4.14) 

fe^oo 

and on the other hand 

lim S (x N+nk ) = lim S (t n {x n „)) = S(r N x), (4.15) 

where the second equality stems from the continuity of the map S and t n . Because 
S is a generalised Lyapunov function, the sequence S (x n ) is convergent. Therefore 
the subsequences (I4.14p and (|4.15|) must have the same limit. We conclude that 
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S(t n x) = S(x) which contradicts Eq. (|4.12l) . Hence x = x*. Since we have shown 
that any convergent subsequence of r n {x) converges to the same limit x*, it follows by 
Lemma l4~3l that T n (x) is converging to x*. Since that holds for arbitrary x, it follows 
that t is mixing. ■ 

Lemma 4.3 Let x n be a sequence in a sequentially compact topological space X such 
that any convergent subsequence converges to x*. Then the sequence converges to a?*. 

Proof We prove by contradiction: assume that the sequence does not converge to x*. 
Then there exists an open neighbourhood 0(x*) of x* such that for all k G N, there is 
a nk such that x nk €" 0(x*). Thus the subsequence x nk is in the closed space X\0(x*), 
which is again sequentially compact. x nk has a convergent subsequence with a limit 
in X\0(x*), in particular this limit is not equal to x*. ■ 

There is an even more general way of defining Lyapunov functions which we state here 
for completeness. It requires the concept of the quotient topology [149]. 

Definition 4.12 Let r : X — ► X be a map on a topological space X. A continuous 
map S : X — > R is called a multi-central Lyapunov function for r around J 7 C X if 
and only if the sequence S(r n (x)) is point-wise convergent for any x G X and if S 
and t fulfil the following three conditions: S is constant on J 7 , r(JT) C JF, and 

S*(x) = lim S{r n (x)) ^ S(x) ^x^T. (4.16) 

n^oo 

For these functions we cannot hope that the orbit is mixing. We can however show 
that the orbit is "converging" to the set T in the following sense: 

Theorem 4.3 (Multi-central Lyapunov function) Let r : X — > X be a continu- 
ous map on a sequentially compact topological space X. Let S : X — > R be a multi- 
central Lyapunov function for r around T. Let (p : X — > X / ' T be the continuous 
mapping into the quotient space (i.e. <f(x) = [x] for x G X\J r and <p(x) = [J-] for 
x G J 7 ). Then f : X jT — ► X jT given by f([x]) = ip (r (</ 7 ~ 1 (N))) is mixing with fixed 
point [J-]. 

Proof First note that f is well defined because ip is invertible on X/F\[F) and 
t(!F) C J 7 , so that f([^ r ]) = [J 7 ]. Since X is sequentially compact, the quotient space 
X / J 7 is also sequentially compact. Note that for O open, f^ 1 (0) = <p (r _1 (v?" 1 (O))) 
is the image of <p of an open set in X and therefore (by definition of the quotient topol- 
ogy) open in X/J 7 . Hence f is continuous. The function SQsc]) : X/J 7 — ► X/J 7 given 
by >S([x]) = S , (9? _1 ([x])) is continuous and easily seen to be a generalised Lyapunov 
function around [J 7 ]. By Theorem 14.21 it follows that f is mixing. ■ 
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4.3.2 Metric spaces 

We now show that for the particular class of compact topological sets which posses a 
metric, the existence of a generalised Lyapunov function is also a necessary condition 
for mixing. 

Theorem 4.4 (Lyapunov criterion) Let r : X — > X be a continuous map on a 
compact metric space X. Then r is mixing with fixed point x* if and only if a generalised 
Lyapunov function around x* exists. 

Proof Firstly, in metric spaces compactness and sequential compactness are equiva- 
lent, so the requirements of Theorem l4.2l are met. Secondly, for any mixing map r with 
fixed point x*, a generalised Lyapunov function around x* is given by S(x) = d(x*,x). 
In fact, it is continuous because of the continuity of the metric and satisfies 

lim S(r n (x)) = d{x*,x*) = < d(x m ,x) = 5(x), (4.17) 

ri— >oo 

where the equality holds if and only x = x*. We call d(x*,x) the trivial generalised 
Lyapunov function. ■ 

Remark 4-1 In the above Theorem we have not used all the properties of the metric. 
In fact a continuous semi-metric (i.e. without the triangle inequality) would suffice. 

The trivial Lyapunov function requires knowledge of the fixed point of the map. There 
is another way of characterising mixing maps as those which bring elements closer to 
each other (rather than closer to the fixed point). 

Definition 4.13 A map r : X — ► X is on a metric space is called a non-expansive 
map if and only if 

d(r(x), r(y)) sC d(x, y) Vx, y e X, (4.18) 
a weak contraction if and only if 

d(r(x), r(y)) < d(x, y) Vx,yeX,x^ y, (4.19) 

and a strict contraction if and only if there exists a k < 1 such that 

d(r(x),r(j/)) < kd(x,y) \/x,y G X . (4.20) 

Remark The notation adopted here is slightly different from the definitions used by 
other Authors [5, 143, 150] who use contraction to indicate our non-expansive maps. 
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Our choice is motivated by the need to clearly distinguish between non-expansive 
transformation and weak contractions. 

We can generalise the above definition in the following way: 

Definition 4.14 A map r : X — ► X on a metric space is called an asymptotic defor- 
mation if and only if the sequence <i(r n (x), r n (y)) converges point-wise for all x, y G X 
and 

lim d(T n (x),T n (y)) ± d(x, y) Vx,yeX,x^ y. (4.21) 

n^oo 

Lemma 4.4 Let r : X — > X be a non-expansive map on a metric space X, and let 

d(T N (x),T N (y)) <d{x,y) Vx,yeX,x^y (4.22) 

for some fixed N £ N. Then r is an asymptotic deformation. Then r is an asymptotic 
deformation. 

Proof The existence of the limit linin^oo d{r n (x) , r n (y)) follows from the monotonic- 
ity and the fact the any metric is lower bounded. ■ 

Remark 4-4 Any weak contraction is an asymptotic deformation (with N = 1). 



Theorem 4.5 (Asymptotic deformations) Let r : X — > X be a continuous 
map on a compact metric space X with at least one fixed point. Then r is mixing 
if and only if r is an asymptotic deformation. 



Proof Firstly assume that r is an asymptotic deformation. Let fixed point 

and define S(x) = d(x*,x). 

lim S(T n (x)) = lim d{x^T n {x)) 

n^oo n— >oo 

= lim d(T n (x*),T n (x)) + d(x*,x) = S(x) Vx/x*, (4.23) 

n— >oo 

hence S(x) is a generalised Lyapunov function. By Theorem 14.21 it follows that r is 
mixing. Secondly, if r is mixing, then 

lim d(T n (x),T n (y)) = d(x*,x*) = / d(x,y) Vx,y£X,x^y, (4.24) 

n— >oo 

so r is an asymptotic deformation. ■ 
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Remark 4-5 Note that the existence of a fixed point is assured if r is a weak contraction 
on a compact space [151], or if the metric space is convex compact [152]. 

As a special case, we get the following result: 

Corollary 4.2 Any weak contraction r on a compact metric space is mixing. 

Proof Since the space is compact r has at least one fixed point. Moreover from 
Lemma l4~4l we know that r is an asymptotic deformation. Then Theorem 03] applies. ■ 

Remark 4-6 This result can be seen as an instance of Banach contraction principle on 
compact spaces. In the second part of the chapter we will present a counterexample 
which shows that weak contractivity is only a sufficient criterion for mixing (see Ex- 
ample H3]). In the context of quantum channels an analogous criterion was suggested 
in [143, 146] which applied to strict contractions. We also note that for weak and 
strict contractions, the trivial generalised Lyapunov function (Theorem 14. 4p is a strict 
Lyapunov function. 

Lemma 14.51 states the ergodic theorem by Birkhoff [153] which, in the context of 
normed vector spaces, shows the equivalence between the definition of ergodicity of 
Eq. (|4.4p and the standard time average definition. 



Lemma 4.5 Let X be a convex and compact subset of a normed vector space, and let 
t : X — > X be a continuous map. If r is ergodic with fixed point x*, then 



1 n 

— Xy(s)=s, . (4.25) 



hm 

n— »oc n + _ 

Proof Define the sequence A n = ^j-j- Yli=o T ^( x )- Let then M be the upper bound 
for the norm of vectors in X ', i.e. M = sup^g^ ||x|| < oo. which exists because X is 
compact. The sequence A n has a convergent subsequence A Uk with limit A. Since r 
is continuous one has lim^oo r(A nk ) = t(A). On the other hand, we have 

r(A ni ) - An. = T nk+1 (x) - x ^ J! v ; " " " < , 4.26 

" V nkJ n fc + l" V ' " " n k + l ^ n k + V K ' 

so the two sequences must have the same limit, i.e. r{A) = A. Since r is ergodic, we 
have A = x* and lim n ^ 00 A n = x* by Lemma I4.31 ■ 

Remark 4-7 Note that if r has a second fixed point y* ^ x*, then for all n one has 
7^1 E"=o r ^(y*) = V*: so E( l- ( ! 4 - 25 P would not apply. 
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4.4 Quantum Channels 

In this Section we discuss the mixing properties of quantum channels [2] which account 
for the most general evolution a quantum system can undergo including measure- 
ments and coupling with external environments. In this context solving the mixing 
problem (14. ip is equivalent to determine if repetitive application of a certain physical 
transformation will drive any input state of the system (i.e. its density matrices) into 
a unique output configuration. The relationship between the different mixing criteria 
one can obtain in this case is summarised in Fig. 14.31 

At a mathematical level quantum channels correspond to linear maps acting on the 
density operators p of the system and satisfying the requirement of being completely 
positive and trace preserving (CPT). For a formal definition of these properties we 
refer the reader to [5, 154, 155]: here we note only that a necessary and sufficient 
condition to being CPT is to allow Kraus decomposition [154] or, equivalently, Stine- 
spring dilation [156]. Our results are applicable if the underlying Hilbert space is 
finite-dimensional. In such regime there is no ambiguity in defining the convergence 
of a sequence since all operator norms are equivalent (i.e. given two norms one can 
construct an upper and a lower bound for the first one by properly scaling the second 
one). Also the set of bounded operators and the set of operators of Hilbert-Schmidt 
class coincide. For the sake of definiteness, however, we will adopt the trace-norm 
which, given the linear operator : TL — > TL, is defined as ||6||i = Tr[V 0t0] with 
Tr[- • • 1 being the trace over TL and 8t being the adjoint of G. This choice is in part 
motivated by the fact [150] that any quantum channel is non-expansive with respect 
to the metric induced^] by || • ||i (the same property does not necessarily apply to other 
operator norms, e.g. the Hilbert-Schmidt norm, also when these are equivalent to 

We start by showing that the mixing criteria discussed in the first half of the chapter 
do apply to the case of quantum channel. Then we will analyse these maps by studying 
their linear extensions in the whole vector space formed by the linear operators of TL. 

4.4.1 Mixing criteria for Quantum Channels 

Let TL be a finite dimensional Hilbert space and let S(TL) be the set of its density 
matrices p. The latter is a convex and compact subset of the larger normed vector 
space C(TL) composed by the linear operators : TL — > TL of TL. From this and from 
the fact that CPT maps are continuous (indeed they are linear) it follows that for 

3 This is just the trace distance d(p,a) — \\p — a\\i. 
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Figure 4.3: Relations between the different properties of a quantum channel. 



a quantum channel there always exists at least one density operator which is a fixed 
point [140]. It also follows that all the results of the previous section apply to quantum 
channels. In particular Lemma [4.11 holds, implying that any mixing quantum channel 
must be ergodic. The following example shows, however, that it is possible to have 
ergodic quantum channels which are not mixing. 

Example Consider the qubit quantum channel r obtained by cascading a com- 
pletely decoherent channel with a NOT gate. Explicitly r is defined by the transfor- 
mations t(|0)(0|) = |1)(1|, t(|1><1|) = |0)(0|, and r(|0)(l|) = r(|l)(0|) = with |0), |1) 
being the computational basis of the qubit. This map is ergodic with fixed point given 
by the completely mixed state (|0)(0| + |l)(l|)/2. However it is trivially not mixing 
since, for instance, repetitive application of r on |0)(0| will oscillate between |0)(0| 
and |1)(1|. 

Theorems 14.51 implies that a quantum channel r : S(Tt) — ► S(Tt) is mixing if and 
only if it is an asymptotic deformation. As already pointed out in the introduction, 
this property is metric independent (as opposed to contractivity) . Alternatively, if the 
fixed point of a quantum channel is known, then one may use the trivial generalised 
Lyapunov function (Theorem 14. 4|) to check if it is mixing. However both criteria 
depend on the metric distance, which usually has no easy physical interpretation. A 
more useful choice is the quantum relative entropy, which is defined as 

H{p, a) = Trp(log p - log a) . (4.27) 

The quantum relative entropy is continuous in finite dimension [157] and can be used 
as a measure of distance (though it is not a metric). It is finite if the support of p is 
contained in the support of a. To ensure that it is a continuous function on a compact 
space, we choose a to be faithful: 

Theorem 4.6 (Relative entropy criterion) A quantum channel with faithful fixed 
point p* is mixing if and only if the quantum relative entropy with respect to p* is a 
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generalised Lyapunov function. 

Proof Because of Theorem 14.21 we only need to prove the second part of the thesis, 
i.e. that mixing channels admit the quantum relative entropy with respect to the fixed 
point, S(p) = H(p,p*), as a generalised Lyapunov function. Firstly notice that the 
quantum relative entropy is monotonic under quantum channels [158,159]. Therefore 
the limit S*(p) = lim n ^oo S (r n (p)) does exist and satisfies the condition S*(p) ^ S(p). 
Suppose now there exists a p such that S*(p) = S(p). Because r is mixing and S is 
continuous we have 

S(p) = S*(p) = lim S(r n (p)) = S(p*) = 0, (4.28) 

n— >oo 

and hence H(p,p*) = 0. Since H(p,a) = if and only if p = a it follows that S is a 
Lyapunov function around p*. ■ 

Another important investigation tool is Corollary 14. 21 weak contractivity of a quantum 
channel is a sufficient condition for mixing. As already mentioned in the previous 
section, unfortunately this not a necessary condition. Here we present an explicit 
counterexample based on a quantum channel introduced in Ref. [140]. 

Example 4-3 Consider a three-level quantum system characterised by the orthogo- 
nal vectors |0),|1),|2) and the quantum channel r defined by the transformations 
r(|2}(2|) = |1)(1|, r(|l)(l|) = r(|0)(0|) = |0)(0|, and r(\i)(j\) = for all i + j. Its 
easy to verify that after just two iterations any input state p will be transformed into 
the vector |0)(0|. Therefore the map is mixing. On the other hand it is explicitly not 
a weak contraction with respect to the trace norm since, for instance, one has 

|| r(|2)(2|) - r(|0)(0|) ||i = || |1><1| - |0)<0| ||i = || |2)(2| - |0)(0| \\t , (4.29) 

where in the last identity we used the invariance of || • ||i with respect to unitary 
transformations. 

4.4.2 Beyond the density matrix operator space: spectral properties 

Exploiting linearity quantum channels can be extended beyond the space S(Tt) of 
density operators to become maps defined on the full vector space C(Ti) of the linear 
operators of the system, in which basic linear algebra results hold. This allows one 
to simplify the analysis even though the mixing property (|4.1|) is still defined with 
respect to the density operators of the system. 
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Mixing conditions for quantum channels can be obtained by considering the struc- 
ture of their eigenvectors in the extended space C{TL). For example, it is easily shown 
that the spectral radius [160] of any quantum channel is equal to unity [140], so its 
eigenvalues are contained in the unit circle. The eigenvalues A on the unit circle (i.e. 
|A| = 1) are referred to as peripheral eigenvalues. Also, as already mentioned, since 
S(Tt) is compact and convex, CPT maps have always at least one fixed point which 
is a density matrix [140]. 

Theorem 4.7 (Spectral gap criterion) Let t be a quantum channel, r is mixing 
if and only if its only peripheral eigenvalue is 1 and this eigenvalue is simple. 

Proof The "if" direction of the proof is a well known result from linear algebra (see 
for example [160, Lemma 8.2.7]). Now let us assume r is mixing towards p*. Let 6 be a 
generic operator in C{TL). Then 6 can be decomposed in a finite set of non-orthogonal 
density operator^, i.e. © = Yli c lPlt with pi £ S(7i) and q complex. Since Tr [p$\ = 1, 
we have have Tr [G] = Yle c i- Moreover since r is mixing we have lmin^oo r n (pi) = p* 
for all £, with convergence with respect to the trace-norm. Because of linearity this 
implies 



If there existed any other eigenvector 6* of r with eigenvalue on the unit circle, then 
lim n ^oo T n (G*) would not satisfy Eq. (|4.30D . ■ 

The speed of convergence can also be estimated by [140] 



where N is the dimensionality of the underlying Hilbert space, k is the second largest 
eigenvalue of r, and Cn is some constant depending only on N and on the chosen 
norm. Hence, for n N the convergence becomes exponentially fast. As mentioned 
in [143], the criterion of Theorem 14.71 is in general difficult to check. This is because 
one has to find all eigenvalues of the quantum channel, which is hard especially in 
the high dimensional case. Also, if one only wants to check if a particular channel 

4 To show that this is possible, consider an arbitrary operator basis of C{7i). If N is the finite 
dimension of Ti the basis will contain N 2 elements. Each element of the basis can then be 
decomposed into two Hermitian operators, which themselves can be written as linear combinations 
of at most JV projectors. Therefore there exists a generating set of at most 27V 3 positive operators, 
which can be normalised such that they are quantum states. There even exists a basis (i.e. a 
minimal generating set) consisting of density operators, but in general it cannot be orthogonalised. 




(4.30) 



r n (p)-p*\\i ^C N n N K 



(4.31) 
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is mixing or not, then the amount of information obtained is much higher than the 
required amount. 

Example 4-4 As an application consider the non mixing CPT map of Example 14.21 
One can verify that apart from the eigenvalue 1 associated with its fixed point (i.e. 
the completely mixed state), it possess another peripheral eigenvalue. This is A = — 1 
which is associated with the Pauli operator |0)(0| — |1)(1|. 



Corollary 4.3 The convergence speed of any mixing quantum channel is exponentially 
fast for sufficiently high values of n. 



Proof From Theorem 14.71 mixing channels have exactly one peripheral eigenvalue, 
which is also simple. Therefore the derivation of Ref. [140] applies and Eq. (|4.31|) 
holds. ■ 



This result should be compared with the case of strictly contractive quantum channels 
whose convergence was shown to be exponentially fast along to whole trajectory [143, 
146]. 



4.4.3 Ergodic channels with pure fixed points 

An interesting class of ergodic quantum channel is formed by those CPT maps whose 
fixed point is a pure density matrix. Among them we find for instance the maps 
employed in the communication protocols discussed in this thesis or those of the pu- 
rification schemes of Refs. [141, 142]. We will now show that within this particular 
class, ergodicity and mixing are indeed equivalent properties. 

We first need the following Lemma, which discusses a useful property of quantum 
channels (see also [161]). 

Lemma 4.6 Let r be a quantum channel and be an eigenvector of r with peripheral 



eigenvalue A = e 1 ^ . Then, given g = Tr V 0t0 > 0, the density matrices p 
and a = V Q^Q/g are fixed points of r. 



Proof Use the left polar decomposition to write = g pU where U is a unitary 
operator. The operator pU is clearly an eigenvector of r with eigenvalue e l<p , i.e. 

T (pU) = XpU . (4.32) 
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Hence introducing a Kraus set {K n } n of r [154] and the spectral decomposition of the 
density matrix p = Pj\i/jj)(ipj\ with pj > being its positive eigenvalues, one gets 

A = Tr[T(pU)tf] = ^pji^K^i^lUKlU^^) , (4.33) 

where the trace has been performed with respect to an orthonormal basis {|^)}^ of 
TL. Taking the absolute values of both terms gives 

|A| = \Y,Pj(^\ K n\^)(^\UKlU^\^)\ 



= ^WU)}VWipJ] = h (4.34) 

where the inequality follows from the Cauchy-Schwartz inequality The last identity 
instead is a consequence of the fact that the transformation f(p) = Ut(U^ pU)U^ is 
CPT and thus trace preserving. Since |A| = 1 it follows that the inequality must be 
replaced by an identity. This happens if and only if there exist e tl3 such that 

^{(falK^j)}* = VPj(^'Kl^) = erty/pjtyjpKlrtlfa) , (4.35) 

for all j,£ and n. Since the \4>i) form a basis of H., and pj > this implies 

{^\Kl = e vd {^ j \UKitf => (^\UKl = e-^ (^\KlU , (4.36) 

for all n and for all the not null eigenvectors of p. This yields 

j n j n 

= e~ m r{p)U (4.37) 



which, replaced in (14.321) gives e~ l ® r(p) = e tip p, whose only solution is e~ tl3 = e tif '. 
Therefore r(p) = p and p is a fixed point of r. The proof for a goes along similar 
lines: simply consider the right polar decomposition of G instead of the left polar 
decomposition. ■ 

Corollary 4.4 Let r be an ergodic quantum channel. It follows that its eigenvectors 
associated with peripheral eigenvalues are normal operators. 
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Proof Let G be an eigenoperator with peripheral eigenvalue e J</3 such that r (0) 



we know that, given g 
'^Q/g must 

a, i.e. eet = ete. 



Tr 



the density matrices 



e tip 0. By Lemma 

p = V 60t/<? and a = V €rf@/g must be fixed points of r. Since the map is ergodic 
we must have p 



Theorem 4.8 (Purely ergodic maps) Let be the pure fixed point of 

an ergodic quantum channel r. It follows that r is mixing. 



Proof We will use the spectral gap criterion showing that |V ; i)(V'i| is the only pe- 
ripheral eigenvector of r. Assume in fact that 6 E L(H.) is a eigenvector of r with 
peripheral eigenvalue, i.e. 

r(e) = e^e. (4.38) 

Prom Lemma l4~6l we know that the density matrix 

p = Veet/s, (4.39) 

with g = Tr V 0^6 > 0, must be a fixed point of r. Since this is an ergodic map we 
must have p = This implies G = g\ipi) (ip2\ , with some normalised vector 

of TL. Replacing it into Eq. (|4.38|) and dividing both terms by g yields r (l^iX^D = 
(-02 1 and 

KViKI^M^DIV^Hi- (4.40) 

Introducing a Kraus set {K n } n of r and employing Cauchy-Schwartz inequality one 
can then write 

1 = K^|r(|^i)(^|)|^2)| = |j;(^il«n|^i><^|^l^2>l (4-41) 



where we used the fact that l^i) is the fixed point of r. Since r is CPT the quantity 
(ip2 1 r ( I ^2) 1) I V'2) is upper bounded by 1. Therefore in the above expression the 
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inequality must be replaced by an identity, i.e. 

(^2|r(|V 2 )(^2|)|^2) = l ^ t(|^ 2 )(^ 2 |) = |V>2>(^2|- (4.42) 

Since r is ergodic, we must have IV^KV^I = l^iKV'il- Therefore oc \ipi){ipi\ which 
shows that is the only eigenvector of r with peripheral eigenvalue of. ■ 

An application of the previous Theorem is obtained as follows. 

Lemma 4.7 Let Mab = Ma ® 1_b + lyl ® Mb be an observable of the composite system 
Ti-A ®7~Lb and T the CPT linear map on TLa of Stinespring form [156] 

t( P ) =Tr B [U(p® |0) B (^|)C/tj ( (4.43) 

(here Trx [• • • ] is the partial trace over the system X , and U is a unitary operator of 
TLa®7~Lb)- Assume that [Mab,U] = and that \4>)b is the eigenvector corresponding 
to a non- degenerate maximal or minimal eigenvalue of Mb- Then r is mixing if and 
only if U has one and only one eigenstate that factorises as \v)a ® \<t>)B- 

Proof Let p be an arbitrary fixed point of r (since r is CPT it has always at least 
one), i.e. Tr B [U (p ® 17+] = p. Since Mab is conserved and Tr^ [M^p] = 

Tica [Mat(p)], the system B must remain in the maximal state, which we have assumed 
to be unique and pure, i.e. 

U(p®\<f>) B (<j>\)rt =p®\<f>)B(<j>\ p7,p®|$ B (0|]=O. (4.44) 

Thus there exists a orthonormal basis {\u k )} k oH1.a®'Hb diagonalising simultaneously 
both U and p ® |0) B (0|. We express the latter in this basis, i.e. p ® \4>)b{4>\ = 
J2kPk\ u k){uk\ with pk > 0, and compute the von Neumann entropy of subsystem B. 
This yields 



= H (\4>) b {<ft\) = H I Tia 



^2pk\uk)(uk\ 




p k H(TTA[\u k )(u k \])(4A5) 



From the convexity of the von Neumann entropy the above inequality leads to a 

contradiction unless Tta [|^fc)(«fc|] = \4>} B (<j)\ for all k. The \u k ) must therefore be 
factorising, 

K) = K)a® \4>)b- (4.46) 

If the factorising eigenstate of U is unique, it must follow that p = |z^)<V| for some 
\u) and that r is ergodic. By Theorem 14.81 it then follows that r is also mixing. If on 
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the other hand there exists more than one factorising eigenstate, than all states of the 
form of Eq. (|4.46|) correspond to a fixed point pt = \vk)(yk\ and r is neither ergodic 
nor mixing. ■ 

Remark 4-3 An application of this Lemma is the protocol for read and write access 
by local control discussed in the next chapter. 

4.5 Conclusion 

In reviewing some known results on the mixing property of continuous maps, we 
obtained a stronger version of the direct Lyapunov method. For compact metric spaces 
(including quantum channels operating over density matrices) it provides a necessary 
and sufficient condition for mixing. Moreover it allows us to prove that asymptotic 
deformations with at least one fixed point must be mixing. 

In the specific context of quantum channels we employed the generalised Lyapunov 
method to analyse the mixing properties. Here we also analysed different mixing 
criteria. In particular we have shown that an ergodic quantum channel with a pure 
fixed point is also mixing. 
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5.1 Introduction 

The unitarity of Quantum Mechanics implies that information is conserved. Whatever 
happens to a quantum system - as long as it is unitary, the original state can in 
principle be recovered by applying the inverse unitary transformation. However it is 
well known that in open quantum systems [121] the reduced dynamics is no longer 
unitary. The reduced dynamics is described by a completely positive, trace preserving 
maps, and we have seen in the last chapter that there are extreme examples, namely 
mixing maps, where all information about the initial state is eventually lost. Where 
has it gone? If the whole system evolves unitary, then this information must have been 
transferred in the correlations between reduced system and environment [162], and/or 
in the environment. We can see that this may be useful for quantum state transfer, in 
particular the case where all information is transferred into the "environment", which 
could be another quantum system (the receiver). A particularly useful case is given by 
mixing maps with pure convergence points, because a pure state cannot be correlated, 
and because we have a simple convergence criterion in this case (Subsection 14.4.31) ■ 
This is an example of homogenisation [138,139]. Furthermore, if the mixing property 
arises from some operations, we can expect that by applying the inverse operations, 
information can also be transferred back to the system. This property was used 
in [137, 163] to generate arbitrary states of a cavity field by sending atoms through 
the cavity. The crucial difference is that in our system control is only assumed to be 
available on a subsystem (such as, for example, the ends of a quantum chain). Hence 
we will show in this chapter how arbitrary quantum states can be written to (i.e. 
prepared on) a large system, and read from it, by local control only. This is similar 
in spirit to universal quantum interfaces [164], but our different approach allows us to 
specify explicit protocols and to give lower bounds for fidelities. We also demonstrate 
how this can be used to significantly improve the quantum communication between two 
parties if the receiver is allowed to store the received signals in a quantum memory 
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before decoding them. In the limit of an infinite memory, the transfer is perfect. 
We prove that this scheme allows the transfer of arbitrary multi-partite states along 
Heisenberg chains of spin-1/2 particles with random coupling strengths. 

Even though the convergence of a mixing map is essentially exponentially fast 
(Corollary 14. 3|) . we still have to deal with infinite limits. Looking at the environment 
this in turn would require to study states on an infinite dimensional Hilbert space, 
and unfortunately this can introduce many mathematical difficulties. We are mainly 
interested in bounds for the finite case: if the protocol stops after finitely many steps, 
what is the fidelity of the reading/writing? Which encoding and decoding operations 
must be applied? By stressing on these questions, we can actually avoid the infinite 
dimensional case, but the price we have to pay is that our considerations become a 
bit technically involved. 



5.2 Protocol 

We consider a tripartite finite dimensional Hilbert space given by TL = TLc®'Hc®'^m- 
We assume that full control (the ability to prepare states and apply unitary transfor- 
mations) is possible on system C and M, but no control is available on system C. 
However, we assume that C and C are coupled by some time-independent Hamilto- 
nian H. We show here that under certain assumptions, if the system CC is initialised 
in some arbitrary state we can transfer ("read") this state into the system M by apply- 
ing some operations between M and C. Likewise, by initialising the system M in the 
correct state, we can prepare ("write") arbitrary states on the system CC. The system 
M functions as a quantum memory and must be at least as large as the system CC. 
As sketched in Fig. 15.11 we can imagine it to be split into sectors Mi, I.e.. 



L 



Hm = (^) T^Mi 



(5.1) 



with 



d\mTLM t = dimHc- 



(5.2) 



For the reading case, we assume that the memory is initialised in the state 



|0)m 



(g)|0>M £ 



(5.3) 



i 
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where |0) can stand for some generic stat^]. Like in the multi rail protocols considered 
in Chapter [3l we let the system evolve for a while, perform an operation, let it evolve 
again and so forth, only that now the operation is not a measurement, but a unitary 
gate. More specifically, at step £ of the protocol we perform a unitary swap Si between 
system C and systems Mg. After the Lth swap operation the protocol stops. The 
protocol for reading is thus represented by the unitary operator 

W = SlUSl^U ■■■SzU---S x U, (5.4) 

where U G C(TC C q) is the time-evolution operator U = exp{—iHt} for some fixed 
time interval t. As we will see in the next section, the reduced evolution of the system 
C under the protocol can be expressed in terms of the CPT map 

T(p o ) = tr c [u(p c ®\0)c(0\)tf] , (5.5) 

where |0)c is the state that is swapped in from the memory. Our main assumption now 
is that r is ergodic with a pure fixed point (which we denote as |0)^) . By Theorem 14.81 
this implies that r is mixing, and therefore asymptotically all information is transferred 
into the memory. 

For writing states on the system, we just make use of the unitarity of W. Roughly 
speaking, we initialise the memory in the state that it would have ended up in after 
applying W if system CC had started in the state we want to initialise. Then we 
apply the inverse of W given by 

jpt = [/t 5l . . . rf Se . . . U ] S L ^U ] S L . (5.6) 

We will see in Section 15.41 how this gives rise to a unitary coding transformation on 
the memory system, such that arbitrary and unknown states can be initialised on the 
system. The reader has probably noticed that the inverse of W is generally unphysical 
in the sense that it requires backward time evolution, i.e. one has to wait negative 
time steps between the swaps. But we will see later how this can be fixed by a simple 
transformation. For the moment, we just assume that is physical. 



1 Later on we will give an example where |0) represents a multi-qubit state with all qubits aligned, 
but here we don't need to assume this. 
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Figure 5.1: The system CC can only be controlled by acting on a (small) subsystem 
C. However system C is coupled to system C by a unitary operator U = exp{—iHt} . 
This coupling can - in some cases - mediate the local control on C to the full system 
CC. In our case, system C is controlled by performing regular swap operations Sg 
between it and a quantum memory Mg. 



5.3 Decomposition equations 



In this section we give a decomposition of the state after applying the protocol which 
will allow us to estimate the fidelities for state transfer in terms of the mixing properties 
of the map r. Let \ip)cc e ^cc be an arbitrary state. We notice that the C component 
of W\ij}) C q\G) m is always \0)c- Therefore we can decompose it as follows 



W\^) C c\0)m = \0)c ® \y/v\0)c\<t>)M + VT^]\A) 



CM 



(5.7) 



with \A)q M being a normalised vector of C and M which satisfies the identity 



?<0|A) 



CM 



. 



(5.8) 
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It is worth stressing that in the above expression n, \4>)m and |A)g M are depending 
on \ip)cc- We decompose acting on the first term of Eq. (|5.7p as 



W j \0) cC \(p)M = \4)cc\0)m + V^f, |A)ccm, 
where \A) C q M is the orthogonal complement of \iP)cc\®)m, i-e. 



Uf(o|A; 



CCM 



. 



(5.9) 



(5.10) 



Multiplying Eq. (|5.9p from the left with ccX^M^] and using the conjugate of Eq. 
(|5.7|) we find that 77 = fj. An expression of 77 in terms of r can be obtained by noticing 
that for any vector |V ; )cc the following identity applies 



c) 



tr c 



U(p c ®\0) c (0\)Ut 



trcM 



with being the reduced density matrix trc YWIcc 
one gets 



US i (\i(>) Cc (Tp\®\0)M(0\)S i tf\ , 

(5.11) 

Reiterating this expression 



trcM 



(5.12) 



with p'^ = trc (IV')cc(^I) ^ ] • Therefore from Eq. (|5.7|l and the orthogonality 
relation (15.81) it follows that 



r, = c (0\r L - 1 (p'c) |0)c, 



(5.13) 



which, since r is mixing, shows that 77 —* 1 for L — ► 00. Moreover we can use Eq. (14.311) 
to claim that 



II 



\c 



(Olr^ 1 (pfe) |0> 



c 



II 



< \\r L - l {p' e )- 10)0(0111! <fl(L-l)«fe 



(5.14) 



where i? is a constant which depends upon d c = dirnT^c and where k g]0, 1[ is the 
second largest eigenvalue of r. 



5.4 Coding transformation 

Here we derive the decoding/encoding transformation that relates states on the mem- 
ory M to the states that are on the system CC. We first apply the above decomposi- 
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tions Eqs. (15.71) and (15.91) to a fixed orthonormal basis {\ipk)cc} °f ^cc> i- e - 



Wl^ed )** = |0)c &> I v^/fe|O)^|0 fc >Af + a/1 - r/ fc |A fc )c, M 
W^O/ccI^a/m = V^ife |^fc)ccl°)M + \A - Vk \ Ak)cCM- 



(5.15) 



Define a linear operator D on TCm which performs the following transformation 



kIM 



(5.16) 



Here \ipk)M are orthonormal vectors of M which represent the states {\ipk)cc} °^^CC 
(formally they are obtained by a partial isometry from CC to M). The vectors \4>k)M 
are defined through Eq. (|5.15l) - typically they will not be orthogonal. We first show 
that for large L they become approximately orthogonal. 



From the unitarity of and from Eq. (|5.15l) we can establish the following identity 

M{4>h\(t>k') M = \/Vk Vk' <W + VVk (1 - Vk') CCM^k^&k^CCM ( 5 - 17 ) 
+ VV' (1 - Vk) CCM(^k\lpk>0)cCM + V( l ~ %)(! - Vk') C c M (^k\^k')cCM ■ 

Defining rjo = min^ r/k it follows for k ^ k' that 

\M((f>k\(f>k')M\ < VVk (1 - Vk') IcCM^kdl^kdcCMl (5.18) 
+ VVk' (1 - Vk) \cCM(^k\lpk'0)cCM\ 
+ ~ Vk)(l ~ Vk') \cCM(^k\^k')cCM\ 

< 2 a/1 - tto + (1 - r/o) < 3a/1 - (5.19) 
Therefore for all /c, /c' the inequality 

\M(<Pk\4>k')M - S k ,k'\ < 3 \Jl - vo (5.20) 

holds. It is worth noticing that, since Eq. (|5.14l) applies for all input states \ip)cci we 
have 

|% - 1| < C (L - l) d c K L ~ l . (5.21) 
Eq. (|5.20|) allows us to make an estimation of the eigenvalues of D^D as 

\\ k -l\<3d cC a/1-770, (5.22) 
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with d, 



cc 



&rmTi C Q. We now take a polar decomposition D = PV of D. V is the 



best unitary approximation to D [160, p 432] and we have 



\D-V\\ 



k 



< 3 d cQ y/l -7/o. 



Therefore 



|£>-F|| 2 < >/3d oC (l-»») 



1/4 



(5.23) 



(5.24) 



which, thanks to Eq. (|5.2ip . shows that D can be approximated arbitrary well by a 
unitary operator V for L — ► oo. 



5.5 Fidelities for reading and writing 

In what follows we will use and F as our reading and writing transformation, 
respectively. In particular, W will be used to recover the input state \ifi)cc °f ^ ne 
chain after we have (partially) transferred it into M through the unitary W (i.e. we 
first act on \ip)cc ® |0)m with W, and then we apply W on M). Vice-versa, in order 
to prepare a state \ip)cc on we ^ rs ^ prepare M into |V>)a/, then we apply to it 
the unitary transformation V and finally we apply W*. We now give bounds on the 
fidelities for both procedures. 

The fidelity for reading the state \4>)m is given by 

F r (V0 = m(^|W Rm V\^) m (5.25) 
where Rm is the state of the memory after W, i.e. 



R M = tr 



cc 



W(\ip)cc{ip\ ® 10)^/(01)^1 = ^ I^a/^I + (1 " »7) °M • (5.26) 



In the above expression we used Eqs. (|5.7|) and (|5.8|) and defined om = tr^[| A)^j^(A|]. 
Therefore by linearity we get 

FM) = V \m(<P\V\iP) m \ 2 + {1~V) m{1>\V* o-m V\^) m > V \ M {<t>\V\i>) M \ 2 ■ (5.27) 
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Notice that 



\m(<P\V\tP) m \ = \m(4>\V -D + D\^) m \ > \mU>\DW>m\ - \m{<P\D - V\^) M \ ■ (5.28) 



Now we use the inequality (|5.24|) to write 

\m(4>\D-V\iI>)m\ < \\D-V\\ 2 <^d cC (1-%) 1/4 • 



(5.29) 



If \iP)m was a basis state \ipk)Mi then \m{<P\D\iP) m\ = 1 by the definition Eq. (I5.16|) 
of D. For generic \ip)M we can use the linearity to find after some algebra that 



y/rj \m{4\D\^)m\ > Vm ~ 3 d cc V 1 ~ Vo ■ 
Therefore Eq. (|5.28p gives 

y/rj \m((I>\V\^)m\ > Vvo -5d cC (l-r/ ) 1/4 . 



By Eq. ([BT271) it follows that 

F r > m - 10 d cC (1 - 7? ) 1/4 . 



(5.30) 



(5.31) 



(5.32) 



The fidelity for writing a state \ip)cc CC is given by 



>cc- 



(5.33) 



A lower bound for this quantity is obtained by replacing the trace over M with the 
expectation value on |0)a/, i.e. 

i^(V0 > cc(i>\M(0\WiV(\i>) M (i>\®\0)cc(U\)V j! W\0)MH)cc 



cc 



M\W) 



cc 



n 



M 



<V|v^|0> 



M 



V \m{4>\V\iIj)m\ 



(5.34) 



where Eqs. (|5.7p and the orthogonality relation (|5.8|) have been employed to derive the 
second identity. Notice that the last term of the inequality (I5.34p coincides with the 
lower bound (|5.27|) of the reading fidelity. Therefore, by applying the same derivation 
of the previous section we can write 



F>rj - 10 d cC (1 - rjo) 



1/4 



(5.35) 
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Figure 5.2: Alice and Bob control the spins Na and Nb interconnected by the spins 
Nr. At time jt Bob performs a swap Sj between his spins and the memory Mj. 



which shows that the reading and writing fidelities converge to 1 in the limit of large 
L. Note that this lower bound can probably be largely improved. 



5.6 Application to spin chain communication 

We now show how the above protocol can be used to improve quantum state transfer 
on a spin chain. The main advantage of using such a memory protocol is that - opposed 
to all other schemes - Alice can send arbitrary multi-qubit states with a single usage 
of the channel. She needs no encoding, all the work is done by Bob. The protocol 
proposed here can be used to improve the performances of any scheme mentioned in 
Section 11.51 and it works for a large class of Hamiltonians, including Heisenberg and 
XY models with arbitrary (also randomly distributed) coupling strengths. 

Consider a chain of spin-1/2 particles described by a Hamiltonian H which conserves 
the number of excitations. The chain is assumed to be divided in three portions 
A (Alice), B (Bob) and R (the remainder of the chain, connecting Alice and Bob) 
containing respectively the first Na spins of the chain, the last Nb spins and the 
intermediate Nr spins, and the total length of the chain is N = Na+Nr+Nr (see Fig 
I5.2p . Bob has access also to a collection of quantum memories M\, • • • , Mj • • • , Ml 
isomorphic with B, i.e. each having dimension equal to the dimension 2 Nb of B. 
We assume that Bob's memory is initialised in the zero excitation state |0)m- Alice 
prepares an arbitrary and unknown state a on her Na qubits. By defining the 
(from Bob's perspective) controlled part of system C = B and the uncontrolled part 
C = AR, we can apply the results of the last sections and get the following 



Theorem 5.1 (Memory swapping) Let H be the Hamiltonian of an open 
nearest-neighbour quantum chain that conserves the number of excitations. If there 
is a time t such that /i,Ar(t) ^ (i.e. the Hamiltonian is capable of transport be- 
tween Alice and Bob) then the state transfer can be made arbitrarily perfect by 
using the memory swapping protocol. 
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PROOF We only have to show that the reduced dynamics on the chain is mixing with 
a pure fixed point. Using the number of excitations as a conserved additive observable, 
we can use the criterion of Lemma l4~7l If there exists exactly one eigenstate \E) of 
factorising form with |0)b, i.e. 



then the reduced dynamics is mixing toward \0)ar- Assume by contradiction that has 
an eigenvector \E)ar / |0) ar which falsifies Eq. (|5.36p . Such an eigenstate can be 
written as 



where a and b are complex coefficients and where the spin just before the section B 
(with position Na + Nr) is in the state |0) for \/j)ar and in the state |1) for \P)ar- 
Since the interaction between this spin and the first spin of section B includes an 
exchange term (otherwise /i,Ar(t)=0 for all t), then the action of H on the second 
term of (|5.37|) yields exactly one state which contains an excitation in the sector 
B. It cannot be compensated by the action of H on the first term of (|5.37j) . But 
by assumption \E)ar ® |0)s is an eigenstate of H, so we conclude that 6 = 0. This 
argument can be repeated for the second last spin of section R, the third last spin, and 
so on, to finally yield \E)ar = |0)ar, as long as all the nearest neighbour interactions 
contain exchange parts. ■ 

Remark 5.1 Theorem 15.11 should be compared to Theorem 13.21 for the multi rail pro- 
tocol. They are indeed very similar. However the current theorem is much stronger, 
since it allows to send arbitrary multi-excitation states, and also to write states back 
onto the chain. It is interesting to note that Lemma 14.71 and Theorem 15.11 indicate a 
connection between the dynamical controllability of a system and its static entangle- 
ment properties. It may be interesting to obtain a quantitative relation between the 
amount of entanglement and the convergence speed. 

Let us now come back to the question raised in Section 15.21 about the operation W' 
being unphysical. As mentioned before, this can be fixed using a simple transfor- 
mation: if the Hamiltonian H fulfils the requirements of Lemma 14.71 then also the 
Hamiltonian —H fulfils them. Now derive the coding transformation V as given in 
Section EU for the Hamiltonian H = —H. In this picture, the reading protocol W is 
unphysical, whereas the writing protocol becomes physical. In the more general case 



3i |A)ar : H\\) A r®\0)b = E\\)ar®\0)b 



(5.36) 



\E) A r ® |0) B = a\n) A R ® |0)b + b\p)AR ® \0)b 



(5.37) 
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where the condition of Lemma 14.71 is not valid, but the map 



r{pc) = tr c U(pc®\0)c(Q\)tf 



(5.38) 



is still ergodic with a pure fixed point, we then require the map 



f{po) = tr c UHpc®\Q)c(Q\)U 



(5.39) 



to be also ergodic with pure fixed point to be able to use this trick. 

5.7 Conclusion 

We have given an explicit protocol for controlling a large permanently coupled system 
by accessing a small subsystem only. In the context of quantum chain communication 
this allows us to make use of the quantum memory of the receiving party to improve 
the fidelity to a value limited only by the size of the memory. We have shown that 
this scheme can be applied to a Heisenberg spin chain. The main advantage of this 
method is that arbitrary multi-excitation states can be transferred. Also, our method 
can be applied to chains that do not conserve the number of excitations in the system, 
as long as the reduced dynamic is ergodic with a pure fixed point. 

It remains an open question how much of our results remain valid if the channel is 
mixing toward a mixed state. In this case, a part of the quantum information will in 
general remain in the correlations between the system and the memory, and it cannot 
be expected that the fidelity converges to one. However, by concentrating only on the 
eigenstate of the fixed point density operator with the largest eigenvalue, it should be 
possible to derive some bounds of the amount of information that can be extracted. 
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6.1 Introduction 



We have mainly discussed two methods for quantum state transfer so far. In the 
first one, multiple chains where used, and in the second one, a single chain was used 
in combination with a large quantum memory. Can we combine the best of the two 
schemes, i.e. is it possible to use only a single chain and a single memory qubit? In this 
chapter we will show that this is indeed the case and that the fidelity can be improved 
easily by applying in certain time-intervals two-qubit gates at the receiving end of the 
chain. These gates act as a valve which takes probability amplitude out of the system 
without ever putting it back. The required sequence is determined a priori by the 
Hamiltonian of the system. Such a protocol is optimal in terms of resources, because 
two-qubit gates at the sending and receiving end are required in order to connect the 
chain to the blocks in all above protocols (though often not mentioned explicitly). At 
the same time, the engineering demands are not higher then for the memory swapping 
protocol. Our scheme has some similarities with [92], but the gates used here are 
much simpler, and arbitrarily high fidelity is guaranteed by a convergence theorem 
for arbitrary coupling strengths and all non-Ising coupling types that conserve the 
number of excitations. Furthermore, we show numerically that our protocol could 
also be realised by a simple switchable interaction. 



6.2 Arbitrarily Perfect State Transfer 

We now show how the receiver can improve the fidelity to an arbitrarily high value 
by applying two-qubit gates between the end of the chain and a "target qubit" of the 
block. We label the qubits of the chain by 1, 2, • • • , N and the target qubit by N + 1 
(see Fig. I6.ip . The coupling of the chain is described by a Hamiltonian H. We assume 
that the Hamiltonian H conserves the number of excitations and that the target qubit 
N + 1 is uncoupled, 

H\N + 1)=0 (6.1) 
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and set the energy of the ground state |0) to zero. For what follows we restrict all 
operators to the N + 2 dimensional Hilbert space 

H = span {\n); n = 0, 1, 2, . . . , N + 1} . (6.2) 

Our final assumption about the Hamiltonian of the system is that there exists a time 
t such that 

f N , t {t) = (N\exp{-itH}\l) + 0. (6.3) 

Physically this means that the Hamiltonian has the capability of transporting from 
the first to the last qubit of the chain. As mentioned in the introduction, the fidelity 
of this transport may be very bad in practice. 



V k 




1 2 N N+l 

Figure 6.1: A quantum chain (qubits 1,2, ■■■ ,N) and a target qubit (N + 1). By 
applying a sequence of two-qubit unitary gates Vk on the last qubit of the chain and 
the target qubit, arbitrarily high fidelity can be achieved. 



We denote the unitary evolution operator for a given time tk as Uk = exp {— it^H} 
and introduce the projector 

P = 1 - |0)(0| - \N)(N\ -\N + 1}{N + 1|. (6.4) 

A crucial ingredient to our protocol is the operator 

V(c,d) = P + |0)(0| +d\N)(N\ +d*\N + 1}(N + 1| 

+c*|JV + l)(JV|-c|JV)(JV + l|, (6.5) 

where c and d are complex normalised amplitudes. It is easy to check that 

VV ] = V*V = 1, (6.6) 

so V is a unitary operator on 7i. V acts as the identity on all but the last two 
qubits, and can hence be realised by a local two-qubit gate on the qubits N and N+l. 
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Furthermore we have VP = P and 

V(c, d) [{c\N) + d\N + 1)}] = | AT + 1). (6.7) 

The operator V(c, d) has the role of moving probability amplitude c from the iVtli 
qubit to target qubit, without moving amplitude back into the system, and can be 
thought of as a valve. Of course as V(c, d) is unitary, there are also states such that 
V(c, d) acting on them would move back probability amplitude into the system, but 
these do not occur in the protocol discussed here. 

Using the time-evolution operator and two-qubit unitary gates on the qubits N and 
N + 1 we will now develop a protocol that transforms the state |1) into \N + 1). Let 
us first look at the action of U\ on |1). Using the projector P we can decompose this 
time-evolved state as 

Ui\l) = PU X \1) + \N){N\U X \1) 

= PC7i|l} + VPl{ci|iV)+di|JV-r-l)}, (6.8) 

where p x = \{N\U\\1)\ , a = (N\Ui\l)/y/pi and di = 0. Let us now consider the 
action of V\ = V(ci,di) on the time-evolved state. By Eq. (|6.7|) it follows that 

VW X |1) = P?7i|l) + v^r|iV + l). (6.9) 

Hence with a probability of pi, the excitation is now in the position N + 1, where it is 
"frozen" (since that qubit is not coupled to the chain. We will now show that at the 
next step, this probability is increased. Applying JJi to Eq. (|6.9p we get 

U 2 V 1 U 1 \1) 

= PU 2 PU l \l) + (N\U 2 PU X \1)\N) + ^\N + 1) 

= PU 2 PUx\l) + v^{c 2 |iV) + d 2 \N + 1)} (6.10) 

with c 2 = (N\U 2 PU X \1) / d 2 = jpiljpi and 

P2 = Pi + |(AT| ^PC/ill) | 2 > Pl . (6.11) 

Applying V 2 = V(c 2 ,d 2 ) we get 

V 2 U 2 V 1 Ui\l) = PU 2 PU 1 \l) + y /pi\N + l). (6.12) 
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Repeating this strategy £ times we get 



,k=i j \k=i j 



(6.13) 



where the products are arranged in the time-ordered way. Using the normalisation of 
the r.h.s. of Eq. (16 . 13|) we get 



Pe = 1 




(6.14) 



From Section 13.51 we know that there exists a t > such that for equal time intervals 
t\ = t2 = ■ ■ ■ = = t we have lim^ooj^ = 1. Therefore the limit of infinite gate 
operations for Eq. (|6.13|) is given by 



lim 



\[v k u k 



IJV + 1) 



(6.15) 



\k=l 



It is also easy to see that lim^oo dg = 1, lim^oo eg = and hence the gates V k 
converge to the identity operator. Furthermore, since V^f7fc|0) = |0) it also follows 



that arbitrary superpositions can be transferred. As discussed in Theorem 14.311 this 
convergence is asymptotically exponentially fast in the number of gate applied (a 
detailed analysis of the relevant scaling can be found in Chapter Equation (I6.15|) 
is a surprising result, which shows that any non-perfect transfer can be made arbitrarily 
perfect by only applying two-qubit gates on one end of the quantum chain. It avoids 
restricting the gate times to specific times (as opposed to the dual rail scheme) while 
requiring no additional memory qubit (as opposed to the memory swapping scheme). 

The sequence that needs to be applied to the end of the chain to perform the 
state transfer only depends on the Hamiltonian of the quantum chain. The relevant 
properties can in principle be determined a priori by preceding measurements and 
tomography on the quantum chain (as discussed in Sect. I2.9| ). 



6.3 Practical Considerations 

Motivated by the above result we now investigate how the above protocol may be 
implemented in practice, well before the realisation of the quantum computing blocks 
from Fig. ll.4l The two-qubit gates V k are essentially rotations in the {|01), |10)} space 
of the qubits N and N + 1. It is therefore to be expected that they can be realised (up 
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to a irrelevant phase) by a switchable Heisenberg or XY type coupling between the 
Nth and the target qubit. However in the above, we have assumed that the gates Vk 
can be applied instantaneously, i.e. in a time-scale much smaller than the time-scale 
of the dynamics of the chain. This corresponds to a switchable coupling that is much 
stronger than the coupling strength of the chain. 
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Figure 6.2: Numerical example for the convergence of the success probability. Sim- 
ulated is a quantum chain of length N = 20 with the Hamiltonian from Eq. (|6.16|) 
(dashed line) and Eq. (|6.17|) with B/J = 20 (solid line). Using the original proto- 
col [1], the same chain would only reach a success probability of 0.63 in the above time 
interval. 



Here, we numerically investigate if a convergence similar to the above results is still 
possible when this assumption is not valid. We do however assume that the switching 
of the interaction is still describable by an instantaneous switching (i.e. the sudden 
approximation is valid). This assumption is mainly made to keep the numerics simple. 
We do not expect qualitative differences when the switching times become finite as 
long as the time-dependent Hamiltonian is still conserving the number of excitations 
in the chain. In fact it has recently been shown that the finite switching time can even 
improve the fidelity [33]. Intuitively, this happens because by gradually decreasing the 
coupling, he not only receives the probability amplitude of the last qubit of the chain, 
but can also "swallow" a bit of the dispersed wave-packed (similar to the situation 
discussed in [92]). 

We have investigated two types of switching. For the first type, the coupling itself 
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is switchable, i.e. 

N-l 

H(t) = JJ2 a n<+i + A(i)o£<r+ +1 + h.c, (6.16) 

n=l 

where A(i) can be or 1. For the second type, the target qubit is permanently coupled 
to the remainder of the chain, but a strong magnetic field on the last qubit can be 
switched, 

N 

H(t) = J ^ <Tn<7n+l + hc - + BA(t)a* N+1 , (6.17) 

n=l 

where again A(t) can be or 1 and B> 1. This suppresses the coupling between the 
iVth and N + 1th qubit due to an energy mismatch. 

In both cases, we first numerically optimise the times for unitary evolution t^ over a 
fixed time interval such that the probability amplitude at the iVth qubit is maximal. 
The algorithm then finds the optimal time interval during which A(t) = 1 such that 
the probability amplitude at the target qubit is increased. In some cases the phases 
are not correct, and switching on the interaction would result in probability amplitude 
floating back into the chain. In this situation, the target qubit is left decoupled and the 
chain is evolved to the next amplitude maximum at the iVth qubit. Surprisingly, even 
when the time-scale of the gates is comparable to the dynamics, near-perfect transfer 
remains possible (Fig 16. 2|) . In the case of the switched magnetic field, the achievable 
fidelity depends on the strength of the applied field. This is because the magnetic field 
does not fully suppress the coupling between the two last qubits. A small amount of 
probability amplitude is lost during each time evolution U^, and when the gain by the 
gate is compensated by this loss, the total success probability no longer increases. 

6.4 Conclusion 

We have seen that by having a simple switchable interaction acting as a valve for 
probability amplitude, arbitrarily perfect state transfer is possible on a single spin 
chain. In fact, by using the inverse protocol, arbitrar)Q states in the first excitation 
sector can also be prepared on the chain. Furthermore, this protocol can easily be 
adopted to arbitrary graphs connecting multiple senders and receivers (as discussed 
for weakly coupled systems in [86]). 

1 Opposed to the method for state preparation developed in the last chapter this allows the creation 
of known states only (as the valve operations Vk depend explicitly on the state that one wants to 
prepare). 
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7.1 Introduction 

An important question that was left open so far is what happens to quantum state 
transfer in the presence of external noise. It is well known from the theory of open 
quantum systems [121] that this can lead to dissipation and decoherence, which also 
means that quantum information is lost. The evolution of a closed quantum system 
is described by the Schrodinger equation 

= -iH\il>). (7.1) 

If a system is very strongly coupled to a environment, the dynamic is completely inco- 
herent and described by some simple rate equations for the occupation probabilities, 

n n 

In the more general case where the dynamic consists of coherent and incoherent parts, 
the evolution can sometimes be expressed as a Lindblad equation [121] 

dtp = Cp (7.3) 

for the reduced density matrix. These three regimes are shown in Fig. 17.11 For 
quantum information theory, coherence is essential [2], and one has to try to isolate 
the quantum chain as much as possible from the environment. In the partially coherent 
regime, typically the quantum behaviour decays exponentially with a rate depending 
on the temperature of the environment. Not surprisingly, this has also been found 
in the context of quantum state transfer [165-167]. From a theoretical point of view 
it is perhaps more interesting to look at the low temperature and strong coupling 
regime, where the dynamics is often non-Markovian [121] and can no longer expressed 
as a simple Lindblad equation. This is also interesting from a practical perspective, 
corresponding to effects of the environment which cannot be avoided by cooling. Here 
we consider a model where the system is coupled to a spin environment through 
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incoherent 




System coupling 



Figure 7.1: Dominant regimes of dynamics depending on the relative strength of the 
system Hamiltonian and the environmental coupling [47]. 



an exchange interaction. This coupling offers the unique opportunity of an analytic 
solution of our problem without any approximations regarding the strength of system- 
environment coupling (in most treatments of the effect of an environment on the 
evolution of a quantum system, the system-environment coupling is assumed to be 
weak) and allows us to include inhomogeneous interactions of the bath spins with the 
system. For such coupling, decoherence is possible for mixed (thermal) initial bath 
states [168,169]. However if the system and bath are both initially cooled to their 
ground states, is there still a non-trivial effect of the environment on the fidelity? In 
this chapter we find that there are two important effects: the spin transfer functions 
(Eq. I1.19|) are slowed down by a factor of two, and destabilised by a modulation of 
|cosGi| , where G is the mean square coupling to the environment. This has both 
positive and negative implications for the use of strongly coupled spin systems as 
quantum communication channels. The spin transfer functions also occur in the charge 
and energy transfer dynamics in molecular systems [47] and in continuous time random 
walks [170] to which our results equally apply. 
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7.2 Model 

We choose to start with a specific spin system, i.e. an open spin chain of arbitrary 
length N, with a Hamiltonian given by 

N-l 

H S = --J2 J £ (Wt+i + Y t Y l+1 ) , (7.4) 

where Jg are some arbitrary couplings and Xg and Yi are the Pauli-X and Y matrices 
for the £th spin. Toward the end of the section we will however show that our results 
hold for any system where the number of excitations is conserved during dynamical 
evolution. In addition to the chain Hamiltonian, each spin £ of the chain interacts with 
an independent bath of Mg environmental spins (see Fig 17. 2p via an inhomogeneous 
Hamiltonian, 

fc=l 




Figure 7.2: A spin chain of length N = 5 coupled to independent baths of spins. 

In the above expression, the Pauli matrices Xg and Yg act on the £th spin of the 

(£) (t) 

chain, whereas X^ and Y" fc act on the kth environmental spin attached to the £th 
spin of the chain. We denote the total interaction Hamiltonian by 

N 

Hi^H®. (7.6) 

The total Hamiltonian is given by H = H$ + Hi, where it is important to note that 
[Hs, Hi] ^ 0. We assume that a homogeneous magnetic field along the z-axis is applied. 
The ground state of the system is then given by the fully polarised state |0,0), with 
all chain and bath spins aligned along the z-axis. The above Hamiltonian describes an 
extremely complex and disordered system with a Hilbert space of dimension 2 N+NM . 
In the context of state transfer however, only the dynamics of the first excitation sector 
is relevant. We proceed by mapping this sector to a much simpler system [171-175]. 
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For £ = 1, 2, . . . , N we define the states 

\£,0) = X e \0,0) (7.7) 

and 

ln./\ = 



M-^E^M (7-8) 

fc=l 



with 



\ fe=i 

It is easily verified that (setting Jq = Jn = 0) 



Me 9 

EM?)- (™> 



fTsl^.O) = -J/_ 1 |£-1,0)-J^+1,0) 

i?s|0,^> = 0, (7.10) 



and 



H^O) = -G £ \0,£) (7.11) 
Hj\0,£) = -G e \£,0). (7.12) 

Hence these states define a 2N— dimensional subspace that is invariant under the 
action of H. This subspace is equivalent to the first excitation sector of a system of 
2N spin 1/2 particles, coupled as it is shown in Fig 17.31 



Figure 7.3: In the first excitation sector, the system can be mapped into an effective 
spin model where the bath spins are replaced by a single effective spin, as indicated 
here for N = 5. 

Our main assumption is that the bath couplings are in effect the same, i.e. Gg = G 
for all £. Note however that the individual number of bath spins Mg and bath couplings 
g k may still depend on £ and k as long as their means square average is the same. 
Also, our analytic solution given in the next paragraph relies on this assumption, but 
numerics show that our main result [Equation (|7.28l) ] remains a good approximation 
if the Gt slightly vary and we take G = (Ge) . Disorder in the vertical couplings is 
treated exactly in the sense that our results hold for any choice of couplings J#. 
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7.3 Results 



In this paragraph, we solve the Schrodinger equation for the model outlined above and 
discuss the spin transfer functions. Firstly, let us denote the orthonormal eigenstates 
of Hs alone by 

H s \4>k} = e k \4>k) (k = l,2...,N) (7.13) 

with 

N 

Mfc)=5>MM). (7.14) 
t=\ 

For what follows, it is not important whether analytic expressions for the eigensystem 
of Hs can be found. Our result holds even for models that are not analytically solvable, 
such as the randomly coupled chains considered in Section ESI We now make an ansatz 
for the eigenstates of the full Hamiltonian, motivated by the fact that the states 

|#) = -L(|40) + (-l) n |0,*>) (n = l,2) (7.15) 

are eigenstates of H\ with the corresponding eigenvalues ±G [this follows directly 
from Eqs. (|7.1ip and (|7.12|) ]. Define the vectors 

N 

with k = 1,2, ... ,N and n = 0,1. The form an orthonormal basis in which we 
express the matrix elements of the Hamiltonian. We can easily see that 

^|^)=-(-lfC7|^) (7.17) 

and 

JV 

H s \*t) = ^^a M \l,0) = f (|*2) + • (7.18) 
Therefore the matrix elements of the full Hamiltonian H = Hs + Hi are given by 

<F6\H\n) = hw (- (-l) n G5 nn , + |) . (7.19) 

The Hamiltonian is not diagonal in the states of Eq. (|7.16p . But H is now block 
diagonal consisting of iV blocks of size 2, which can be easily diagonalised analytically. 
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The orthonormal eigenstates of the Hamiltonian are given by 

\Ek) = ell {((-If A fc - 2G) + £fc |*l)} (7.20) 

with the eigenvalues 

E%= l -{e k + {-l) n A k ) (7.21) 

and the normalisation 

cun = ^((-l) n A k -2G) 2 + el (7.22) 

where 

A k = ^4G 2 + e 2. (7.23) 

Note that the ansatz of Eq. (|7.16p that put H in block diagonal form did not depend 
on the details of H$ and H\ . The methods presented here can be applied to a much 
larger class of systems, including the generalised spin star systems (which include an 
interaction within the bath) discussed in [175]. 



After solving the Schrodinger equation, let us now turn to quantum state transfer. 
The relevant quantity [1,92] is given by the transfer function 

/jv,i(t) = <J\r,0|exp{-ifft}|l,0) 

= ^ex P {-iE^}(E£|l,0>(iV,0|if>. 

k.n 

The modulus of /jv,i(£) is between (no transfer) and 1 (perfect transfer) and fully 
determines the fidelity of state transfer. Since 



</,(«) = c-„ 1 {((-l)"A fc -2G')^0|* Q fc ) + e,(A0|*^} 

C kn 

V2 



%((-!)" A fc -2G + e k )a u 



we get 



In At) = (7.24) 
2 it (("I)" A, - 2G) 2 + el 



2 

a k\0*kN • 



Eq. (|7.24p is the main result of this section, fully determining the transfer of quantum 
information and entanglement in the presence of the environments. In the limit G — > 0, 
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we have ~ e k and /jv,i(i) approaches the usual result without an environment, 



(7.25) 



In fact, a series expansion of Eq. (17. 24ft yields that the first modification of the transfer 
function is of the order of G 2 , 



G 2 ^2 a k ia* kN 



, , 1 it\ 1 

exp{-ite k \ — 2 + ~2 

' 4 £kj 4 



(7.26) 



Hence the effect is small for very weakly coupled baths. However, as the chains get 
longer, the lowest lying energy ei usually approaches zero, so the changes become more 
significant (scaling as 1/efc). For intermediate G, we evaluated Eq. (|7.24|) numerically 
and found that the first peak of the transfer function generally becomes slightly lower, 
and gets shifted to higher times (Figures 17.41 and I7.5p . A numeric search in the 
coupling space { Jg, £ = 1, . . . , N — 1} however also revealed some rare examples where 
an environment can also slightly improve the peak of the transfer function (Fig l7.6|) . 




Figure 7.4: The absolute value of the transport function /jv,l(i) of an uniform spin 
chain (i.e. Ji = 1) with length N = 10 for three different values of the bath coupling 
G. The filled grey curve is the envelope of the limiting function for G 3> £fc/2 given by 
|/°(|)|. We can see that Eq. (I7.28P becomes a good approximation already at G = 4. 

In the strong coupling regime G>£fc/2,we can approximate Eq. (|7.23|) by w 2G. 
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Time [1/J] 

Figure 7.5: The same as Fig. I7.41 but now for an engineered spin chain [i.e. Jg = 
yJ{N — £)] as in Subsection II. 5.1[ For comparison, we have rescaled the couplings 
such that Yli Jz 1S the same as in the uniform coupling case. 



Inserting it in Eq. (|7.24p then becomes 



In At) 



-iGt 



exp <j -ite k - } a kl a kN + 



co S (Gt)f° N A 



a kN 



(7.27) 



This surprisingly simple result consists of the normal transfer function, slowed down 
by a factor of 1/2, and modulated by a quickly oscillating term (Figures EH and [73]). 
We call this effect destabilisation. Our derivation actually did not depend on the 
indexes of f(t) and we get for the transfer from the nth to the mth spin of the chain 
that 

(7.28) 



fn,m{t) « COs(Gt)/ (- 



It may look surprising that the matrix f n ^ m is no longer unitary. This is because 
we are considering the dynamics of the chain only, which is an open quantum sys- 
tem [121]. A heuristic interpretation of Eq. (|7.28|) is that the excitation oscillates back 
and forth between the chain and the bath (hence the modulation), and spends half 
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of the time trapped in the bath (hence the slowing). If the time of the maximum of 
the transfer function \f® m (t)\ for G = is a multiple of 7r/2G then this maximum is 
also reached in the presence of the bath. We remark that this behaviour is strongly 
non-Markovian [121]. 

Finally, we want to stress that Eq. (|7.28j) is universal for any spin Hamiltonian that 
conserves the number of excitations, i.e. with [Hs,^2e Ze] = 0- Thus our restriction 
to chain-like topology and exchange couplings for H$ is not necessary. In fact the only 
difference in the whole derivation of Eq. (|7.28p for a more general Hamiltonian is that 
Eq. (|7.10|) is replaced by 

H S \£,0) = ^VK,0). (7.29) 

v 

The Hamiltonian can still be formally diagonalised in the first excitation sector as in 
Eq. (|7.14p . and the states of Eq. (|7.2QI) will still diagonalise the total Hamiltonian 
Hs + Hj. Also, rather than considering an exchange Hamiltonian for the interaction 
with the bath, we could have considered a Heisenberg interaction [176], but only for 
the special case where all bath couplings g k are all the same [177]. Up to some 
irrelevant phases, this leads to the same results as for the exchange interaction. 
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Figure 7.6: A weakly coupled bath may even improve the transfer function for some 
specific choices of the Jg. This plot shows the transfer function |/jv,i(i)| for N = 10. 
The couplings Jg were found numerically. 
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7.4 Conclusion 



We found a surprisingly simple and universal scaling law for the spin transfer functions 
in the presence of spin environments. In the context of quantum state transfer this 
result is double-edged: on one hand, it shows that even for very strongly coupled 
baths quantum state transfer is possible, with the same fidelity and only reasonable 
slowing. On the other hand, it also shows that the fidelity as a function of time 
becomes destabilised with a quickly oscillating modulation factor. In practice, this 
factor will restrict the time-scale in which one has to be able to read the state from 
the system. The results here are very specific to the simple bath model and do not hold 
in more general models (such as these discussed in [165, 167], where true decoherence 
and dissipation takes place). What we intended to demonstrate is that even though 
a bath coupling need not introduce decoherence or dissipation to the system, it can 
cause other dynamical processes that can be problematic for quantum information 
processing. Because the effects observed here cannot be avoided by cooling the bath, 
they may become relevant in some systems as a low temperature limit. 
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Our research on quantum state transfer with spin chains has taken us on a journey 
from a very practical motivation to quite fundamental issues and back again. On one 
hand, our results are quite abstract and fundamental, and have related state transfer 
to number theory, topology and quantum convergence. On the other hand, we have 
developed schemes which are simple and practical, taking into account experimental 
hurdles such as disorder and restricted control. While the multi rail scheme and 
the memory swapping scheme will probably become useful only after much further 
progress in experimental QIT, the dual rail scheme and in particular the valve scheme 
have some good chances to be realised in the near future. 

State transfer with quantum chains has become an area of large interest, with more 
than seventy articles on the subject over the last three years. The most important 
goal now is an experiment that demonstrates coherent transfer on a short chain (say 
of length N > 5). Such an experiment is not only useful building a quantum com- 
puter, but also from a fundamental perspective. For instance, the violation of a Bell- 
inequality between distant entangled solid state qubits would be a milestone in the 
field. Since this requires a very high transfer fidelity, the design of such an experiment 
would probably require system dependent theoretical research on how to overcome 
specific types of noise and how to improve the fidelity for specific Hamiltonians. 
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1.1 In areas of universal control, quantum states can easily be transferred 
by sequences of unitary swap gates Sj^ between nearest neighbours. . 

1.2 Schematic layout of a quantum computer. The solid arrows represent 
the flow of quantum information, and the dashed arrows the flow of 
classical information 

1.3 Permanently coupled quantum chains can transfer quantum states with- 
out control along the line. Note that the ends still need to be control- 
lable to initialise and read out quantum states 

1.4 Small blocks (grey) of qubits (white circles) connected by quantum 
chains. Each block consists of (say) 13 qubits, 4 of which are connected 
to outgoing quantum chains (the thick black lines denote their nearest- 
neighbour couplings). The blocks are connected to the macroscopic 
world through classical wires (thin black lines with black circles at their 
ends) through which arbitrary unitary operations can be triggered on 
the block qubits. The quantum chains require no external control. . . 

1.5 A quantum chain consisting of N = 20 flux qubits [34] (picture and 
experiment by Floor Paauw, TU Delft). The chain is connected to four 
larger SQUIDS for readout and gating 

1.6 Minimal fidelity p(t) for a Heisenberg chain of length N = 50 

1.7 Snapshots of the time evolution of a Heisenberg chain with N = 50. 
Shown is the distribution \f n ,i(t)\ 2 of the wave- function in space at 
different times if initially localised at the first qubit 

1.8 Mean and variance of the state |1) as a function of time. Shown is 
the case N = 50 with the y-axis giving the value relative to the mean 
N/2 + 1 and variance (N 2 — 1)/12 of an equal distribution ^2 \ n )- 

1.9 Approximation of the transfer amplitude for N = 50 around the first 
maximum by Bessel and Airy functions [1,61] 
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1.10 pm(T) as a function of T for different chain lengths. The solid curve is 
given by 1.82(2T) -2 / 3 and corresponds to the first peak of the proba- 
bility amplitude (Eq. [L29D \24 

1.11 Quantum capacity, entanglement of formation (EOF), a lower bound 
for the entanglement of distillation (EOD) and the averaged fidelity as 
a function of p(t). We also show the corresponding chain length which 
reaches this value as a first peak and the classical threshold 3 — 2\^2. 
The explicit expression for the quantum capacity plotted here is given 
in [54], and the lower bound of the entanglement of distillation will be 
derived in Section lcT4l [28 

1.12 Snapshots of the time evolution of a quantum chain with engineered 
couplings (11.471) for N = 50. Shown is the distribution of the wave- 
function in space at different times if initially localised at the first qubit 
(compare Fig. CCU) [30 

2.1 Two quantum chains interconnecting A and B. Control of the systems 

is only possible at the two qubits of either end [36 

2.2 Quantum circuit representation of conclusive and arbitrarily perfect 
state transfer. The first gate at Alice's qubits represents a NOT gate 
applied to the second qubit controlled by the first qubit being zero. The 
qubit \iPa)i on the left hand side represents an arbitrary input state at 
Alice's site, and the qubit \tPb)i represents the same state, successfully 
transferred to Bob's site. The i^-gate represents the unitary evolution 

of the spin chains for a time interval of ti [37 

2.3 Semilogarithmic plot of the joint probability of failure P(£) as a function 
of the number of measurements £. Shown are Heisenberg spin- 1/2- 
chains with different lengths N. The times between measurements ti 
have been optimised numerically [39 

2.4 Time t needed to transfer a state with a given joint probability of failure 
P across a chain of length N. The points denote exact numerical data, 

and the fit is given by Eq. (I2.15P [41 

2.5 The minimal joint probability of failure P(£) for chains with length iV 
in the presence of amplitude damping. The parameter J/Y of the curves 

is the coupling of the chain (in Kelvin) divided by the decay rate (ns~ l ). [44 

2.6 Two disordered quantum chains interconnecting A and B. Control of 

the systems is only possible at the two qubits of either end [46 
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2.7 The absolute values of the transition amplitudes /j\r,i(t) and <?jv,i(i) 
for two Heisenberg chains of length N = 10. The couplings strengths 
of both chains were chosen randomly from the interval [0.8J, 1.2J] . 
The circles show times where Bob can perform measurements without 
gaining information on a and (5 [47 

2.8 The relevant properties for conclusive transfer can be determined by 
measuring the response of the two systems at their ends only [51 

2.9 Time t needed to transfer a state with a given joint probability of failure 
P across a chain of length N with uncorrelated fluctuations of A = 0.05. 
The points denote numerical data averaged over 100 realisations, and 
the fit is given by Eq. (|2.53|) . This figure should be compared with Fig. 

El where A = M 

2.10 Most general setting for conclusive transfer: A black box with two inputs 
and two outputs, acting as an amplitude damping channel defined by 
Eqs. (I2~54l and (1235D \55 

2.11 A simple counterexample for a vertically symmetric system where dual 

rail encoding is not possible. The black lines represent exchange couplings. [56 

2.12 An example for a vertically symmetric system where dual rail encoding 

is possible. The black lines represent exchange couplings of equal strength. [57 

3.1 Schematic of the system: Alice and Bob operate M chains, each con- 
taining A" spins. The spins belonging to the same chain interact through 
the Hamiltonian H which accounts for the transmission of the signal in 
the system. Spins of different chains do not interact. Alice encodes the 
information in the first spins of the chains by applying unitary trans- 
formations to her qubits. Bob recovers the message in the last spins of 

the chains by performing joint measurements [59 

3.2 Example of our notation for M = 5 chains of length N = 6 with K = 2 
excitations. The state above, given by |0}i(g)|3}2<8>|0)3(g)|l}4<g)|0)5, has 
excitations in the chains m\ = 2 and mi = 4 at the horizontal position 
m = 3 and ri2 = 1. It is in the Hilbert space 7i(Se) corresponding 
to the subset Sq = {2, 4} (assuming that the sets Si are ordered in 
a canonical way, i.e. Si = {1,2}, 52 = {1,3} and so on) and will 
be written as |(3, 1);6)). There are („) = 10 different sets Sg and the 
number of qubits one can transfer using these states is log 2 10 « 3. The 
efficiency is thus given by R « 3/5 which is already bigger than in the 
dual rail scheme 61 
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3.3 Optimal rates (maximisation of Eq. (13.451 with respect to M) for the 
Multi Rail protocol. Shown are three curves corresponding to different 
values of the joint probability of failure P one plans to achieve [71 

4.1 Schematic examples of the orbits of a ergodic and a mixing map. ... [72 

4.2 Relations between topological spaces [149]. The space of density matri- 
ces on which quantum channels are defined, is a compact and convex 
subset of a normed vectors space (the space of linear operators of the 
system) which, in the above graphical representation fits within the set 

of compact metric spaces [76 

4.3 Relations between the different properties of a quantum channel. ... [85 

5.1 The system CC can only be controlled by acting on a (small) subsystem 
C. However system C is coupled to system C by a unitary operator 
U = exp {—iHt} . This coupling can - in some cases - mediate the local 
control on C to the full system CC. In our case, system C is controlled 
by performing regular swap operations Si between it and a quantum 
memory Mi [96 

5.2 Alice and Bob control the spins Na and Nb interconnected by the 
spins Nr. At time jt Bob performs a swap Sj between his spins and 

the memory Mj llOl 

6.1 A quantum chain (qubits 1, 2, • • • ,N) and a target qubit (N + 1). By 
applying a sequence of two-qubit unitary gates Vk on the last qubit of 

the chain and the target qubit, arbitrarily high fidelity can be achieved. Il05 

6.2 Numerical example for the convergence of the success probability. Sim- 
ulated is a quantum chain of length N = 20 with the Hamiltonian from 
Eq. ([S36D (dashed line) and Eq. (liUTD with B/J = 20 (solid line). Us- 
ing the original protocol [1], the same chain would only reach a success 
probability of 0.63 in the above time interval Il08 

7.1 Dominant regimes of dynamics depending on the relative strength of 

the system Hamiltonian and the environmental coupling [47] Illl 

7.2 A spin chain of length N = 5 coupled to independent baths of spins. . Il 12 

7.3 In the first excitation sector, the system can be mapped into an effective 
spin model where the bath spins are replaced by a single effective spin, 

as indicated here for N = 5 113 



124 



List of Figures 



7.4 The absolute value of the transport function /jv,i(*) of an uniform spin 
chain (i.e. Ji = 1) with length N = 10 for three different values of the 
bath coupling G. The filled grey curve is the envelope of the limiting 
function for G > e k /2 given by |/°(f)|. We can see that Eq. flPED 
becomes a good approximation already at G = 4 Ill6 

7.5 The same as Fig. 17.41 but now for an engineered spin chain [i.e. = 
y/£(N — £)] as in Subsection 11.5.11 For comparison, we have rescaled 
the couplings such that Ji is the same as in the uniform coupling 
case 1 11 7 

7.6 A weakly coupled bath may even improve the transfer function for some 
specific choices of the Ji. This plot shows the transfer function |/jv,i(i)| 

for N = 10. The couplings Ji were found numerically Il 18 
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2.1 The total time t and the number of measurements M needed to achieve 
a probability of success of 99% for different fluctuation strengths A 
(uncorrelated case). Given is the statistical mean and the standard 
deviation. The length of the chain is N = 20 and the number of random 
samples is 10. For strong fluctuations A = 0.1, we also found particular 
samples where the success probability could not be achieved within the 
time range searched by the algorithm 

2.2 The total time t and the number of measurements M needed to achieve 
a probability of success of 99% for different correlations c between the 
couplings [see Eq. (|2.50|) and Eq. (|2.5ip ]. Given is the statistical mean 
and the standard deviation for a fluctuation strength of A = 0.05. The 
length of the chain is N = 20 and the number of random samples is 20. 
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